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ABSTRACT 



We derive an empirical effective temperature and bolomctric luminosity calibra- 
tion for G and K dwarfs, by applying our own implementation of the InfraRed Flux 
Method to multi-band photometry. Our study is based on 104 stars for which we have 
excellent BV {RI)c J H K s photometry, excellent parallaxes and good metallicities. 
' Colours computed from the most recent synthetic libraries (ATLAS9 and MARCS) 

^SJ ' are found to be in good agreement with the empirical colours in the optical bands, 

but some discrepancies still remain in the infrared. Synthetic and empirical bolomctric 
corrections also show fair agreement. 
\ A careful comparison to temperatures, luminosities and angular diameters ob- 

. tained with other methods in literature shows that systematic effects still exist in the 

calibrations at the level of a few percent. Our InfraRed Flux Method temperature 
scale is 100 K hotter than recent analogous determinations in the literature, but is 
in agreement with spectroscopically calibrated temperature scales and fits well the 
\^ ' colours of the Sun. Our angular diameters are typically 3% smaller when compared to 

, other (indirect) determinations of angular diameter for such stars, but are consistent 

with the limb-darkening corrected predictions of the latest 3D model atmospheres and 
also with the results of asteroseismology. 

Very tight empirical relations are derived for bolomctric luminosity, effective tem- 
O ' perature and angular diameter from photometric indices. 

\ We find that much of the discrepancy with other temperature scales and the 

C/3 . uncertainties in the infrared synthetic colours arise from the uncertainties in the use 

' of Vega as the flux calibrator. Angular diameter measurements for a well chosen set 

of G and K dwarfs would go a long way to addressing this problem. 

•i-H . 

. Key words: stars : fundamental parameters - atmospheres - late-type ~ Hertzsprung- 

^-^ ' Russell (HR) diagram - techniques : photometry - infrared : stars. 



1 INTRODUCTION 

Temperatures, luminosities and radii are amongst the basic 
physical data against which models of stellar structure and 
evolution are tested. In this paper, we address one particu- 
lar area of the Hertzsprung-Russell (HR) diagram, that of 
the G and K dwarfs — focusing particularly on stars for 
which the effects of stellar evolution have been negligible 
or nearly negligible during their lifetimes. Somewhat higher 
mass stars than those considered here have been extensively 
studied historically because their luminosities may be used 
to infer stellar ages. Lower mass stars of late-G to K spec- 
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tral types, have been neglected to some extent, probably 
because there have been few secondary benefits in getting 
stellar models right for these stars, and the lack of good 
parallax, diameter and other data. This situation is rapidly 
changing: firstly, Hipparcos has provided the requisite par- 
allax data; secondly, interferometric techniques are making 
the measurement of diameters for such small stars likely 
to be routine within a few years; thirdly, we have our own 
motivation in the form of an ongoing project to follow the 
chemical history of the Milky Way from lower mass stars, 
for which we can infer indirectly the amount of their con- 
stituent helium via stellar luminosity (Jimenez et al 2003). 
In order to achieve our goals, we need accurate and homoge- 
neous effective temperatures and luminosities for our G+K 
dwarfs. 



© 0000 RAS 



2 Casagrande, Portinari and Flynn 



The effective temperature of a stellar surface is a mea- 
sure of the total energy, integrated over all wavelengths, ra- 
diated from a unit of surface area. Since its value is fixed by 
the luminosity and radius, it is readily calculated for theo- 
retical stellar models, and as one of the coordinates of the 
physical HR diagram, it plays a central role in discussions of 
stellar evolution. Most observations, however, provide spec- 
troscopic or photometric indicators of temperature that are 
only indirectly related to the effective temperature. 

The effective temperature of lower main sequence stars 
is not easy to determine and different measurement tech- 
niques are still far from satisfactory concordance (e.g. 
Mishenina & Kovtyukh 2001, Kovtyukh et al. 2003). At 
present, spectroscopic temperature determinations return 
values that are some 100 K hotter that most of the other 
techniques. Even in a restricted and thoroughly studied re- 
gion like that of the Solar Analogs, effective temperature 
determinations for the same star still differ significantly, by 
as much as 150 K (Soubiran & Triaud, 2004). Models pre- 
dict values that are some 100 K hotter than those measured 
(Lebreton et al. 1999 and references therein). 

In this work, we apply the Infrared Flux Method to 
multi-band BV{RI)cJHKs photometry of a carefully se- 
lected sample of G and K dwarf stars. We compare observed 
with synthetic broad-band colours computed from up-to- 
date (ID) model atmospheres. Such models are then used to 
estimate the missing flux needed to recover bolometric lumi- 
nosities from our photometry. Other than the high quality 
of the observational data, the strength of this work relies 
on using very few basic assumptions: these are the adopted 
Vega absolute calibration and zero-points. This also makes 
clearer the evaluation of possible errors and/or biases in the 
results. Both the absolute calibration and the zero-points 
are expected to be well known and the latest generations of 
model atmospheres produce realistic fluxes for a wide range 
of temperatures, gravities and abundances (Bessell 2005) so 
that the adopted model and calibration are at the best level 
currently available. 

The paper is organized as follows: in Section 2 we de- 
scribe our sample; we compare different libraries of model 
atmospheres with observational data in Section 3 and in Sec- 
tion 4 we present our implementation of the Infrared Flux 
Method along with the resulting temperature scale. In Sec- 
tion 5 we test our scale with empirical data for the Sun and 
solar analogs and with recent interferometric measurements 
of dwarf stars. The comparison with other temperature de- 
terminations is done in Section 6 and in Section 7 and 8 we 
give bolometric calibrations and useful tight relations be- 
tween angular diameters and photometric indices. We con- 
clude in Section 9. 



2 SAMPLE SELECTION AND OBSERVATIONS 

The bulk of the targets has been selected from the sample of 
nearby stars in the northern hemisphere provided by Gray 
et al. (2003). Our initial sample consisted of 186 G and K 
dwarfs with Hipparcos parallaxes for which the relative error 
is less than 14%. Accurate BV{RI)c broadband photometry 
has been obtained for the bulk of the stars at La Palma, 
with additional photometry adopted from Bessell (1990a), 
Reid et al. (2001) and Percival, Salaris & Kilkenny (2003). 



All our stars have accurate spectroscopic metallicities and 
for most of them also a-element abundances. Our sample is 
neither volume nor magnitude limited, but it gives a good 
coverage of the properties of local G and K dwarfs from low 
to high metallicity (Figure 0. The sample has been cleaned 
of a number of contaminants, reducing it to a final sample 
of 104 stars, as described in the following sections. 

2.1 Removing double/multiple and variable stars 

Particular attention has been paid to removing unresolved 
double/multiple stars. These stars primarily affect spectro- 
scopic and photometric measurements, making the system 
appear brighter and redder. The Hipparcos catalogue was 
used to make a prior exclusion of 'certain', 'possible' and 
'suspected' (i.e. good, fair, poor, uncertain and suspected) 
multiple systems, on the basis of the quality of the Mult- 
Flag field in the Hipparcos catalogue. 

Some unrecognized double/multiple stars almost cer- 
tainly remain in the sample, since a few stars are found in 
a 'binary main sequence' above the bulk of the main se- 
quence, as Kotoneva, Flynn & Jimenez (2002) suspected in 
such data but were unable to prove. A new technique for 
identifying binaries/multiples has been introduced by Wie- 
len et al. (1999), termed the method, which indicates the 
presence of multiple unresolved systems by using the differ- 
ence between the near-instantaneous direction of the proper 
motion of the star(s) as measured by the Hipparcos satel- 
lite, and the direction of the ground based proper motion 
measured over much longer time scales. Most of the sus- 
pect stars in the Kotoneva, Flynn & Jimenez (2002) sample 
turned out to be A/x multiples, and we have excluded those 
stars which are likely A/x multiples also from the present 
sample (A/i classifications were kindly made for our sample 
stars by Christian Dettbarn in Heidelberg). This reduced 
the sample to 134 stars. 

We have retained widely separated binary systems in 
the sample where components could be studied separately. 
These were retained by making a prior cross-checking of our 
sample with the catalogue of wide-binary and multiple sys- 
tems of Poveda et al. (1994) and then verifying the sepa- 
ration of such systems by inspecting our own or SIMBAD 
images. 13 such stars have been retained. 

Even if dwarfs stars are not expected to show signs of 
strong variability, the existence of ~ 100 or more Hippar- 
cos measurements per star spread over several years, to- 
gether with the excellent temporal stability of the magnitude 
scales, allows the detection of variability at the level of few 
hundredths of a magnitude (Brown et al. 1997) in Hippar- 
cos targets. We use the Hipparcos classification to remove 
so-called 'duplicity induced variables', 'periodic variables', 
'unsolved variables' and 'possible micro- variables'. This re- 
moved a total of 26 stars (however 8 of these were already 
removed as non-single stars). 

The final number of stars satisfying the above require- 
ments was 129. 



2.2 Broad-Band Photometric Observations 

To recover accurate bolometric fluxes and temperatures 
for the stars, we have obtained accurate and homogeneous 
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Figure 1. Distribution of \og{g) and [Fe/H] for the 104 dwarf stars in the basic sample, collected from the literature. As noted in Section 
lOand BX^ adopt ion of an average log(g) of 4.5 for all the stars is sufficiently accurate for our purposes. All stars but one (HD 108564) 
with [Fe/H] < — 1 have direct a-element measurements. 



Johnson- Cousins BV{RI)c and JHKs photometry for all 
the 186 stars in our initial sample. 



2.2.1 Johnson- Cousins photometry 

For most of the stars in our sample with declination north 
of 5 = —25°, we have made our own photometric observa- 
tions from April to December 2004. Observations were done 
from Finland in full remote mode, using the 35-cm telescope 
piggybacked on the Swedish 60-cm telescope located at La 
Palma in the Canary Islands. An SBIG charge-coupled de- 
vice was used through all the observations. Johnson- Cousins 
BV{RI)c colours were obtained for all stars. 

Standard stars were selected from Landolt (1992) and 
among E-regions from Cousins (1980). Although E-regions 
provide an extremely accurate set of stars, there are some 
systematic differences between the Landolt and the SAAO 
system (Bessell, 1995). Therefore, we have used the Landolt 
(1992) standards placed onto the SAAO system by Bessell 
(1995). 

Ten to twenty standard stars were observed each night, 
bracketing our program stars in colour and airmass. When- 
ever possible program stars were observed when passing 
the meridian, in order to minimize extrapolation to zero- 
airmass. Only if the standard deviation between our cal- 
ibrated and the tabulated values for the standards was 
smaller than 0.015 mag in each band was the night was 
considered photometric and observations useful. In addition, 
only program stars for which the final B — V scatter (ob- 
tained averaging five frames) was smaller than 0.015 mag 
were considered usefully accurate. Scatter in the other bands 
was usually smaller. We expect our photometry to have ac- 
curacies of 0.010-0.015 mag on average. 

In addition to our observations, we have gathered 
BV{RI)c photometry of equivalent or better precision from 
the literature : Bessell (1990a) (who also include measure- 
ments from Cousins (1980)), Reid et al. (2001) and Perci- 
val, Salaris & Kilkenny (2003). For stars in common with 
these authors, we have found excellent agreement between 
the photometry, with a scatter of the order of 0.01 mag in 
all bands, and zero-point shifts between authors of less than 
0.01 mag. This is more than adequately accurate for our 
study. 



2. 2. 2 2MASS photometry 

Infrared JHKs photometry for the sample has been taken 
from the 2MASS catalogue. The uncertainty for each ob- 
served 2MASS magnitude is denoted in the catalogue by 
the flags : flags "h_" and "k_msigcom", and is the 

complete error that incorporates the results of processing 
photometry, internal errors and calibration errors. Some 
of our stars are very bright and have very high errors in 
2MASS. We use 2MASS photometry only if the sum of the 
photometric errors is less than 0.10 mag (i.e. "j_"-|-"h_" + 
"k_msigcom"< 0.10). The flnal number of stars is thus re- 
duced from 129 to 104. For our flnal sample the errors in J 
and Ks bands are similar, with a mean value of 0.02 mag, 
whereas a slightly higher mean error is found in H band 
(0.03 mag). 

2.3 Abundances 

We have gathered detailed chemical abundances for our sam- 
ple stars from the wealth of on-going surveys, dedicated to 
investigate the chemical composition of our local environ- 
ment as well as to the host stars of extra-solar planets (see 
references in Table 1). The internal accuracy of such data is 
usually excellent, with uncertainties in the order of 0.10 dex 
or less. However, abundance analysis for late-type dwarfs can 
still be troublesome in some cases (e.g. AUende Prieto et al. 
2004). We are aware that, in gathering spectroscopy from 
different authors, the underlying temperature scale used to 
derive abundances can differ by as much as 50-150 K, which 
would translate into a [Fe/H] error of ~ 0.1 dex (e.g. As- 
plund 2003; Kovtyukh et al. 2003). Therefore, a more con- 
servative error estimate for our abundances would be ~ 0.15 
dex. However, as we show in Section r4.4.2l the InfraRed Flux 
Method, which we employ to recover the fundamental stellar 
parameters, is only weakly sensitive to the adopted metal- 
licity and uncertainties of ±0.15 dex do not bear heavily on 
our main results. 

The best measured elemental abundance in our dwarf 
stars is usually iron (i.e. [Fe/H]) whereas for theoretical mod- 
els the main metallicity parameter is the total heavy-element 
mass fraction, [M/H]. For most of the stars in our sample, 
the spectra provide measurement not only for [Fe/H], but 
also for the a-elements, which dominate the global metallic- 
ity budget. For stars with a-element estimates, we compute 
[M/H] (Yi et al., 2001): 
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Landscape table to go here. 

Table 1. 
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[M/H] = [Fe/H] + log(0.694/c. + 0.306) (1) 
where 

= lOl^/F^l (2) 

is the enhancement factor and [a/Fe] has been com- 
puted by averaging the a-elements. The older formula by 
Salaris et al. (1993) give similar results with a difference 
smaller than 0.02 dex in [M/H] even for the most a-enhanced 
stars ([a/Fe] - 0.4). 

For stars for which the a elements were not available 
we have estimated their contribution from the mean locus 
of the [a/Fe] vs. [Fe/H] relation from the analytical model 
of Pagel & Tautvaisiene (1995). There were 34 such stars in 
our final sample compared to 70 for which a estimates were 
directly available. 



2.4 Reddening corrections 

Interstellar absorption and reddening must be taken into 
account for a correct derivation of stellar parameters, but 
these effects are negligible for our sample stars, as we discuss 
in this section. 

The distance distribution of the sample peaks around 
~ 30 pc, the most distant object being located at a dis- 
tance of ~ 70 pc. For distances closer than 50 pc the po- 
larimetric approach is extremely sensitive and can be used 
as a lower limit even to the expectedly small amounts of 
dust (at least for anisotropic particles). Tinbergen (1982) 
and Leroy (1993b) have confirmed the complete depletion 
of dust within ~ 40 pc from the Sun. Using the catalogue of 
Leroy (1993a) we have found polarimetric measurements for 
21 stars out of the 129 single and non-variable stars selected 
in Secti0n r2.ll The mean percentage polarization is 31 x 10~^ 
which using the Serkowski, Mathewson & Ford (1975) con- 
version factor corresponds to a E{B — V) of 0.0034 mag. 

The Stromgen H/3 index offers an alternative method 
for assessing the reddening of individual stars in our sam- 
ple, especially for the most distant ones; 59 stars in our 
bona fide single and non-variable sample have H/3 measure- 
ments (Hauck & Mermilliod 1998) from which E{b — y) was 
derived using the intrinsic-colour calibration of Schuster & 
Nissen (1989) plus a small offset correction as noted by Nis- 
sen (1994). The reddening distribution for these 59 stars 
peaks around E{b — y) = 0.004, retaining the same mean 
value (and never exceeding 0.015 for a single star) when we 
restrict the sample to distances further than 50 pc. E{B — V) 
can then be calculated using the standard extinction law or 
the relation E{B — V) — 1.35E{b — y) derived by Crawford 
(1975). Once the reddening E{B — V) is known, the extinc- 
tion at any given wavelength can be determined using the 
standard extinction law. Given that the standard error in 
the Schuster & Nissen calibration is on the order of 0.01 
in E{b — y) and observational errors in uvby measurements 
are possible, the implied reddening corrections are below the 
noise level. In addition to this there are indications (Knude 
1979; Vergely et al. 1998) that interstellar reddening is pri- 
marily caused by small dust clouds causing it to vary in steps 
of 0.01-0.03 mag which confirm that the small corrections 
found are consistent with the assumption of zero reddening 
for the sample as a whole. 



3 EMPIRICAL COLOURS VERSUS MODEL 
COLOURS 

Theoretical stellar models predict relations between physi- 
cal quantities such as effective temperature, luminosity and 
stellar radius. From the observational point of view these 
quantities are, in general, not directly measurable and must 
be deduced from observed broad-band colours and magni- 
tudes. In this context, empirical colour-temperature rela- 
tions (Bessell 1979; Ridgway et al. 1980; Saxner & Ham- 
marback 1985; Di Benedetto & Rabbia 1987; Blackwell & 
Lynas-Gray 1994, 1998; Alonso et al. 1996a; Rami'rez & 
Melendez 2005a) and colour-bolometric correction relations 
(Bessell & Wood 1984; Malagnini et al. 1986; Bell & Gustafs- 
son 1989; Blackwell & Petford 1991a; Alonso et al. 1995) are 
normally used. Such calibrations were usually restricted to a 
specific type and/or population of stars, the only exception 
being Alonso et al. (1995, 1996a) who give empirical rela- 
tions for F0-K5 dwarfs with both solar and sub-solar metal- 
licity. Recently Ramirez & Melendez (2005a) have improved 
and extended the Alonso et al. (1996a) colour-temperature- 
metallicity relations. 

In any case, even empirical calibrations make use of 
model atmospheres to some extent, though the model depen- 
dence is expected to be small (see Section r4.3|l . For example, 
the work of Alonso et al. (1996b) was based on Kurucz' 
(1993) spectra, whereas the recent Ramirez & Melendez 
(2005b) calibration is based on the original Kurucz spectra 
(i.e. without empirical modifications) taken from Lejeune et 
al. (1997). 

In this section we compare a number of stellar models 
(ATLAS9, MARCS, BaSel 3.1 and BaSel 2.1.) to our em- 
pirical data, mainly in two colour planes. In general, colour- 
colour plots of the stars against all the model sets look quite 
satisfactory. 

3.1 ATLAS9, MARCS and BaSeL spectral 
libraries 

We have tested synthetic spectra from the ATLAS9- 
ODFNEW (CasteUi & Kurucz 2003), MARCS (Gustafsson 
et al. 2002), BaSel 2.1 and 3.1 (Lejeune et al. 1997 and West- 
era et al. 2002, respectively) libraries. All grids are given in 
steps of 250 K. Notice that, since we are working with dwarf 
stars, we assume log{g) — 4.5 throughout. This assump- 
tion is reasonable given that determinations of log{g) have 
a typical uncertainty of 0.2 dex or more either by requiring 
Fel and Fell lines to give the same iron abundance or by 
using Hipparcos parallaxes (eg. Bai et al. 2004). This un- 
certainty covers the expected range in \og{g) on the main 
sequence (see Figure 0. In any case, a change of ±0.5 dex 
in the surface gravity implies differences that never exceed a 
few degrees in derived effective temperature, as will be seen 
( Section HX^ . 

3.1.1 ATLAS9-0DFNEW models 

The ATLAS9-ODFNEW models calculated by CasteUi & 
Kurucz include improvements in the input physics, the up- 
date of the solar abundances from Anders & Grevesse (1989) 
to Grevesse & Sauval (1998) and the inclusion of new molec- 
ular line lists for TiO and H2O. The metallicities cover 
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Figure 2. Optical and infrared colours for the ATLAS9-ODFNEW from Castelli & Kurucz (2003) compared to empirical colours for G 
and K dwarfs. The colours are shown as a function of {V — Kg) for [M/H] equal to +0.5 (black line), +0.0 (red line), —0.5 (cyan line), 
— 1.0 (green line), —1.5 (blue line), —2.0 (yellow line). Points correspond to observed colours for the sample stars in the range [M/H] 
> 0.25 (black), -0.25 < [M/H] ^ 0.25 (red), -0.75 < [M/H] ^ -0.25 (cyan), -1.25 < [M/H] ^ -0.75 (green), -1.75 < [M/H] -1.25 
(blue). The metallicities given for the model are solar-scaled, whereas [M/H] for the stars has been computed using eq. Q. A typical 
error bar for the points is shown in the upper left of each plot. The model and empirical colours are generally in very good agreement. 



— 2.5 ^ [M/H] ^ +0.5 with solar-scaled abundance ratios. A 
microturbulent velocity ^ = 2 km s~^ and a mixing length 
parameter of 1.25 are adopted. The extension of these mod- 
els to include also a-enhanced chemical mixtures is ongoing. 
However, we do not expect large differences among differ- 
ent models as a result of a-enhancement, as explained in 
Section 14.31 In the remainder of the paper we refer to the 
ATLAS9-0DFNEW model simply as ATLAS9. 



3.1.2 MARCS models 

The new generation of MARCS models (Gustafsson et al. 
1975; Plez et al. 1992) now includes a much-improved treat- 
ment of molecular opacity (Plez 2003; Gustafsson et al. 
2002), chemical equilibrium for all neutral and singly ionised 
species as well as some doubly ionised species, along with 
about 600 molecular species. The atomic lines are based on 
VALD (Kupka et al. 1999). The models span the metal- 
licity range —5.0 ^ [Fe/H] ^ +1.0. We have used models 



with solar relative abundances for [Fe/H] ^ and enhanced 
a-elements abundances for [Fe/H] < 0. The fraction of a- 
elements is given by the model and [M/H] has been com- 
puted using eq. Q. 



3.1.3 The BaS el libraries 

The BaSel libraries are based on a posteriori empirical cor- 
rections to hybrid libraries of spectra, so as to reduce the 
errors of the derived synthetic photometry. BaSel 2.1 (Leje- 
une et al. 1997) is calibrated using solar metallicity data 
only and it is known to be less accurate at low metallici- 
ties ([Fe/H] < 1) especially in the ultraviolet and infrared 
(Westera et al. 2002). The colour calibration has recently 
been extended to non-solar metallicities by Westera et al. 
(2002). Surprisingly, they found that a library that matches 
empirical colour-temperature relations does not reproduce 
Galactic globular-cluster colour-magnitude diagrams, as a 
result of which they propose two different versions of the 
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Figure 3. As Figure|2]but for MARCS models, although the metallicity range is different, since MARCS model are 'o-enhanced'. Lines 
of constant [M/H] correspond to +0.5 (black line), +0.0 (red line), —0.35 (cyan line), —0.69 (green line), —1.19 (blue line), —1.69 (yellow 
line). Points correspond to observed colours for the sample stars in the range [M/H] > 0.25 (black), —0.175 < [M/H] ^ 0.25 (red), 
-0.52 < [M/H] < -0.175 (cyan), -0.94 < [M/H] < -0.52 (green), -1.44 < [M/H] < -0.94 (blue) and [M/H] ^ -1.44 (yellow). Notice 
that MARCS model are a— enhanced for metallicities below the solar and [M/H] for stars and models have been computed using eq. Q. 



library. In what follows we have used the library built to 
match the empirical colour-temperature relations. For both 
BaSel libraries, it is obvious that the transformations used 
do not correct the physical cause of the discrepancies. How- 
ever, it is worthwhile to check whether they do give better 
agreement with the empirical relations. 

All libraries have been tested in the range 4250 ^ 
Tc// < 6500 K and -2.0 ^ [M/H] +0.5, the intervals 
covered by our sample of stars. 



3.2 Comparison of model and empirical colours 

All our colour-colour comparisons are as a function of 
{V — Ks)- This approach highlights the behaviour of colours 
as function of temperature and also make more straightfor- 
ward the comparison with the empirical colour-temperature 
relations in the Section r3.3l {V — Ks) has a large range com- 
pared to its observational uncertainty and it is known to be 
metallicity insensitive (e.g. Bell & Gustafsson 1989 and Fig- 
ure [TSll. 



All models are found to be in good overall agreement 
with the observations. In particular, the improvement in the 
input physics of the latest ATLAS9 and MARCS models 
show excellent agreement with the optical data (Figures |5| 
and and as good as the semi-empirical libraries in the 
other bands, so that we do not find any specific reason to 
prefer the use of semi-empirical libraries. 

We expect the IR colours to be very important both 
in observational terms and in theoretical modelling. Re- 
cently, Fremaux et al. (2006) have presented a detailed com- 
parison of theoretical spectra (from the NeMo library) in 
H band with observed stellar spectra at resolution around 
3000. They show that in the infrared range, although the 
overall shape of the observed flux distribution is matched 
reasonably well, individual features are reproduced by the 
theoretical spectra only for stars earlier than mid F type. For 
later spectral types the differences increase and theoretical 
spectra of K type stars have systematically weaker line fea- 
tures than those found in observations. They conclude that 
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these discrepancies stem from incomplete data on neutral 
atomic (and to a minor extent to molecular) lines. 

The synthetic optical colours agree well with observed 
ones and models also predict the correct metallicity depen- 
dence. In the infrared, models fit satisfactorily the overall 
trend given by the observed colours, however the scatter in 
the data prevent firm conclusions. \n J — H and J — Ks 
the model colours seem to be slightly offset by ~ 0.05 mag, 
whereas in H — Ks the disagreement is somewhat smaller 
though it depends on the models considered and it seems to 
increase going to cooler temperatures. 

3.3 Empirical colour-temperature relations 
compared to observed and synthetic 
photometry 

We now compare empirical colour-temperature calibrations 
with synthetic and observed photometry. To this purpose, 
we use the calibration recently proposed by Ramirez & 
Melendez (2005a). 

In Fig. |S] we plot the difference between the tempera- 
ture of the model atmospheres (BaSel 3.1 and ATLAS9) and 
the temperature expected for their colours according to the 



empirical Ramirez & Melendez relations. Interestingly, all 
models fail to set onto the empirical temperature scale be- 
ing at least 100 K hotter in optical and infrared colours. The 
disagreement is particularly strong at low metallicity as no- 
ticed by Ramirez & Melendez. The semi-empirical libraries 
of Lejeune et al. (1997) and Westera et al. (2002), though 
in principle optimized to reproduce colour-temperature re- 
lations, provide only a marginal improvement with respect 
to the latest purely theoretical ATLAS9 and MARCS li- 
braries. The temperature offset can be translated into a shift 
in colours needed to set the models onto the empirical scale. 
This is done for the ATLAS9 model in the last row of Figure 
|S] in terms of colours, the mismatch corresponds to ~ 0.1 
mag. 

Empirical calibrations depend on metallicities and 
colour indices only. The adopted absolute calibration can 
not be the main cause of the disagreement since all colours 
(in the optical) scale accordingly. In the optical the adopted 
Vega model that sets the absolute calibration has been ac- 
curately tested by Bohlin & Gilliland (2004) (see Appendix 
A). Although we adopt a different absolute calibration in 
the optical and in the infrared (see Appendix A), this could 
at most affect the colour-metallicity-temperature relations 
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Figure 5. Effective temperatures recovered by means of the Ramirez & Melendez empirical calibration are compared with the effective 
temperature of the corresponding BaSel 3.1 (first row) and ATLAS9 (second row) models used to generate synthetic photometry. C ~ ? 
indicate the colours B — V , V — H and V — Kg. The notation AT((^ — ^) means temperature found by means of the {C, — calibration 
minus the effective temperature of the model atmosphere used to generate the synthetic C, and 5 colours. The following metallicities are 
plotted: -1-0.0 (red), —0.5 (cyan), —1.0 (green), —1.5 (blue) and —2.0 (yellow). The dashed area is the standard deviation of the empirical 
calibration in the corresponding colours. The third row shows the shift in colours required to set the ATLAS9 models on the empirical 
relation. 



in V ~ J, V — H and V — Ks- Even using in the infrared 
the same absolutely calibrated Vega model adopted in the 
optical does not eliminate the disagreement and offsets of 
100 K still persist. 

We have also verified that the choice of zero-points plays 
a negligible role: differences in the optical zero-points based 
on Vega or on the average of Vega and Sirius (see Appendix 
A) lead to mean temperature differences of 10 to 30 K. Be- 
sides being much smaller than the uncertainty in the colour- 
temperature relations, such differences go in the direction of 
worsening the disagreement. 

As to the IR zero-points, we checked their effects by 
generating IR synthetic magnitudes in the TCS and Johnson 
system (where the Vega zero-points are different from those 
deduced from 2MASS photometry) and compared the model 
predictions to the Alonso et al. (1996a) empirical calibration 
in these bands; typical discrepancies of more than 100 K still 
persist. 

The causes for the discrepancy must therefore be model 
deficiencies and/or inaccuracies in the adopted empirical re- 
lations. 

In fact we caution that also empirical relations may hide 
some inadequacies. For instance, the Ramirez & Melendez 
calibration predicts slightly bluer colours for the Sun than 
the recent completely empirical determination of Holmberg, 



Flynn & Portinari (2006). We can also combine and invert 
the Ramirez & Melendez' colour-temperature calibrations, 
to derive the corresponding empirical colour-colour relations 
and compare them to the observed colours for our sample 
of stars. Interestingly the agreement is not always good: the 
solar metallicity is slightly offset with respect to the data 
in {R — I)c band and in the infrared tends to oscillate (see 
Figure |5J . This underlines that there is still room for im- 
provement also in the empirical relations. 



4 AN IMPLEMENTATION OF THE IRFM 

In this section we have used BV{RI)cJHKs photometry 
as the basic observational information to derive bolomet- 
ric ffuxes and effective temperatures of our 104 G and K 
dwarfs. Our approach follows the InfraRed Flux Method 
(IRFM) used in the extensive work of Alonso, Arribas & 
Martinez-Roger (1995, 1996b). Ours is a new independent 
implementation of the IRFM, applied to our 104 G and K 
type dwarfs. 

Our implementation differs from Alonso et al (1995, 
1996), as we base its absolute calibration on a synthetic 
spectrum of Vega (as described in detail in Appendix A) 
rather than on semi-empirical measurements (Alonso, Ar- 
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Figure 6. Lines represent the colour-colour relations derived from the empirical colour-temperature calibrations of Ramirez &; Melendez 
(2005a). Points correspond to our sample of stars. Colours and metallicity ranges for lines and points are same as in Figurel^ except that 
metallicities of -1-0.5 have been excluded since the empirical relations do not hold in that range. In the optical the empirical relations 
fit the observed colours with an accuracy comparable to the theoretical ones (compare with Figure |2J. However in the infrared, for 
metallicities around the solar, empirical relations oscillate in reproducing the observed trend. This is especially interesting in light of the 
puzzling infrared colours of the Sun and the shifts required to cool down the temperature scale as discussed in Section |5] 



ribas & Martinez- Roger, 1994). The use of absolutely cali- 
brated synthetic spectra rather than ground-based measure- 
ments traces back to Blackwell et al. (1990, 1991a), who 
concluded that atmospheric models of Vega offered higher 
precision than observationally determined absolute calibra- 
tions. This approach has been adopted in the extensive work 
of Cohen and collaborators and also Bessell, Castelli & Plez 
(1998), who concluded that model atmospheres are more re- 
liable than the near-infrared absolute calibration measure- 
ments. Models are nowadays sophisticated enough to war- 
rant detailed comparisons between observation and theory 
and the latest space based measurements confirm the valid- 
ity of the adopted absolute calibration (see Appendix A). 
We also differ from the recent work of Ramirez & Melendez 
(2005b) in using the most recent model atmospheres and 
absolute calibration available. Furthermore we use only di- 
rect observational data, explicitly avoiding the use of colour 
calibrations from previous studies to infer the bolometric 
luminosity. 

Following Alonso et al. (1995), the flux outside the 
wavelength range covered by our photometry has been es- 
timated using model atmospheres. Since the percentage of 
F^,ol measured in the directly observed bands range from 
~ 70% to ~ 85% depending on the star, the dependence of 
the estimated bolometric flux on the model is small. 

In Section 14.11 and 14.21 we give a detailed description 
of our implementation of the IRFM, mainly following the 
formulation and terminology adopted by Alonso et al. (1995, 
1996b). 

Our procedure can be summarized as follows. Given the 



metallicity, the surface gravity (assuming \og{g) = 4.5 for 
all our dwarfs, see Section and an initial estimate for 
the temperature of a star, we have interpolated over the 
grid of model atmospheres to find the spectrum that best 
matches these parameters. This spectrum is used to esti- 
mate that fraction of the bolometric fiux outside our fil- 
ters {BV{RI)cJHKs), i.e. the 'bolometric correction'. The 
bolometric flux is determined from the observations, includ- 
ing the bolometric correction. A new effective temperature 
Tcff can be computed by means of the IRFM. This temper- 
ature is used for a second interpolation over the grid, and 
the procedure is iterated until the temperature converges to 
within 1 K (typically within 10 iterations). 

We have tested the results using ATLAS9, MARCS, 
BaSel 2.1 and 3.1 libraries. A detailed discussion on the the 
dependence on the adopted model is given in Section [4.31 

In the following sections we describe our implementa- 
tion of the IRFM effective temperature scale. 

4.1 Bolometric fluxes 

As mentioned above, we use grids of synthetic spectra to 
bootstrap our IRFM. The grids of synthetic spectra all have 
a resolution of 250 K in temperature whereas the resolution 
in metallicity depends on the library, but with typical steps 
of 0.25 or 0.50 dex. 

For any given star of overall metallicity [M/H] (c.f. Eq. 
Q, we interpolate over the grid of synthetic spectra in the 
following way. First, we use the To// : {V — K) calibration 
of Alonso et al. (1996a) to obtain an initial estimate of the 
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effective temperature T^fffi of the star. Then we Unearly 
interpolate in log(Te//) bracketing our temperature estimate 
Tefffl at two fixed values of metallicity, which bracket the 
measured [M/H] of the star. A third linear interpolation 
is finally done in metallicity in order to obtain the desired 
synthetic spectrum. 

Having obtained the spectrum F{\) that better matches 
the physical parameters of the star, for any given band 
(running from B to Ks), we convolve it through the trans- 
mission curve rj(A) of the filter and associate the resulting 
J^^ (model) with its effective wavelength Ac/ / (see Appendix 
B). 

We then compute the fiux covered by the passbands B 
to Ks for the model star {Tb-Ks (model)). When the latter 
is compared with the bolometric flux for the same model star 
(.?^Soi (model)) we can obtain an estimate of the fraction of 
the flux C encompassed by our fllters: 

^ ^ Fb-Ks (model) 
TboI (model) 

We now have the correction factor 1/C for the missing 
flux of a given star, we use its observed BV(RI)cJHKs 
magnitudes (m^) to calculate the fiux as it arrives on the 
Earth: 

.Fc (Earth) = jrf* (Earth) 10-"'*<™«-'"c"\ (4) 

where .T^^*'' (Earth) is the absolute calibrated flux on the 
Earth of the standard star and mj*'* is its observed magni- 
tude. The observed magnitudes and the absolute calibration 
of the standard star are those given in Table EH and El re- 
spectively and play a key role in determining both bolomet- 
ric flux and effective temperature as we discuss in Section 

The flux at the Earth for each band is once again as- 
sociated with the corresponding effective wavelength of the 
star, and a simple integration leads to J-B-Xg (Earth). The 
latter is then divided by the C factor as deflned in eq. @ in 
order to obtain the bolometric flux measured on the Earth 
.7^s„, (Earth). 

4.2 Effective Temperatures 

The effective temperature Tc// of a star satisfies the Stefan- 
Boltzman law TboI = cT^ff, where TboI is the bolometric 
flux on the surface of the star. Only the bolometric flux on 
the Earth is measurable, and we must take into account the 
angular diameter {9) of the star. Thus: 

J^so, (Earth) = GT^ff. (5) 

The way to break the degeneracy between effective tem- 
perature and angular diameter in the previous equation is 
provided by the InfraRed Flux Method (Blackwell & Shal- 
lis 1977; Blackwell, Shalhs & Selby 1979; Blackwell, Petford 
& Shallis 1980). The underlying idea relies on the fact that 
whereas the bolometric flux depends on both angular di- 
ameter and effective temperature (to the fourth power), the 
monochromatic flux in the infrared at Earth - Tx^^ (Earth) 
- depends on the angular diameter but only weakly (roughly 
to the flrst power) on the effective temperature: 

.^A:K(Earth)= 0(re//,(7,A:H), (6) 



where (fi{Tcff,g, Am) is the monochromatic surface flux 
of the star. The ratio J-boI (Earth) /jFajpj (Earth) defines what 
is known as the observational _R- factor (Roba), where the 
dependence on 8 is eliminated. In this sense the IRFM can 
be regarded as an extreme example of a colour method for 
determining the temperature. By means of synthetic spectra 
it is possible to define a theoretical counterpart (Rtheo) on 
the surface of the star, obtained as the quotient between the 
integrated fiux {aT^jj) and the monochromatic fiux in the 
infrared .T^Am (model). 

The basic equation of the IRFM thus reads: 

TboI (Earth) _ cr T^f f 
T\^^ (Earth) J^Xi^ (model) 

and can be immediately rearranged to give the effective 
temperature, eliminating the dependence on 6. 

The monochromatic fiux (Earth) is obtained from 
the infrared photometry using the following relation: 

.^Am (Earth) = g(Am) .^Ij^ (Earth) lO-O-*^'"™-™™ (8) 

where ^(Air) is a correction factor to determine 
monochromatic flux from broad-band photometry, mm 
is the infrared magnitude of the target star, mm'* and 
•^Am (Earth) are the magnitude and the absolute monochro- 
matic flux of the standard star in the same infrared band (see 
Table IXTl and IX2l . For the standard and the target star Xm 
refers to their respective effective wavelength in the given 
infrared band, as we discuss in more detail, together with 
the definition of the g-factor, in Appendix B. 

The IRFM is often applied using more than one infrared 
wavelength and different T^ff are obtained for each wave- 
length. Ideally, the derived temperatures should of course 
be independent of the monochromatic wavelength used. In 
our case we have used J, H, Ks photometry and the corre- 
sponding effective wavelengths. 

We have tested the systematic differences in tempera- 
tures, luminosities and diameters obtained when only one 
infrared band at a time is used for the convergence of the 
IRFM. The results are shown in Figure |7| 

The J band returns temperatures that are systemati- 
cally cooler than those found from Ks , which translates into 
greater angular diameters. The H band, returns tempera- 
tures that are systematically hotter and therefore smaller 
angular diameters than in Ks- These systematic differences 
do not appear to depend on the temperature range, at least 
between 4500-6500 K. 

The systematic offsets could be traced back to the ab- 
solute calibration in different bands. We further point out 
that for the sake of this test, the convergence on T^ff in 
the IRFM is required in one band only and in Figure |7| the 
differences are thus exaggerated. The discrepancy can reach 
a few percent in temperatures and diameters, whereas for 
bolometric luminosities is usually within 2% though it in- 
creases with increasing temperature (Figure |7J. This is due 
to the fact that the derived bolometric luminosity is con- 
strained from the full optical and infrared photometry and 
absolute calibration, whereas for temperature and angular 
diameter, when relying on one IR band only, the absolute 
calibration in that band strongly affects the results. 

Alonso et al. (1996b) concluded that the consistency of 
the three infrared bands they used (J, H, K) was good over 
4000 K, but decided to use only H and K below 5000 K 
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Figure 7. Comparison of effective temperatures, angular diam- 
eters and bolometric luminosities obtained in the three infrared 
bands used. Angular diameters have been scaled using the stellar 
parallaxes from Hipparcos, although the parallaxes do not affect 
the comparison since all bands scale accordingly. The luminosi- 
ties of the stars are in good agreement across all the IR bands 
(because luminosity is also tied down by the optical data as dis- 
cussed in the text) . Offsets of this size are likely to stem from the 
difficulties in setting the zero points of the absolute calibration. 
Horizontal and diagonal lines with slope 1 are intended to guide 
the eye. 



and only K below 4000 K and basically the same is done by 
Ramirez & Melendez (2005b). We have decided to use the 
temperature returned in all three bands also below 5000 K, 
since systematic differences are not too large even below 
5000 K, the cooler temperatures in J band being compen- 
sated by the hotter ones in H band. 

At each iteration the temperature used for the conver- 
gence is the average of the three IR values weighted with 
the inverse of their errors "h_," and "kjnsigcom". This 
temperature is then used to select a new model interpolat- 
ing over the grid of synthetic spectra as described in Section 
14.11 and a new bolometric flux and temperature are then 
computed. The procedure is iterated until the average tem- 
perature converges within 1 K. As can be appreciated from 
Figure |5]the systematic differences between the three bands 
are much reduced. Ideally the ratio of the temperatures de- 
termined in the three bands should be unity. The three mean 
ratios T{Ks)/T{J), T{Ks)/T{H) and T{H)/T{J) together 
with their standard deviations, are 1.0069 {a = 0.77%), 
0.9943 (cr = 0.82%) and 1.0126 {a = 0.89%) and this also 
confirms the quality of the adopted absolute calibration. The 
mean standard deviation for the three infrared temperatures 
is 39 K, reflecting the different band sensitivity to the tem- 
perature, the effect of the absolute calibration in different 
bands as well as the observational errors in the infrared pho- 
tometry. 
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Figure 8. Differences as a function of temperature between the 
final adopted values and those obtained in each ( = J, H, Kg band 
when the convergence is done averaging the three IR values. 



4.3 Internal accuracy and dependence on the 
adopted library 

We have probed the accuracy and the internal consistency of 
the procedure described in Sections 14.11 and 14.21 generating 
synthetic magnitudes in the range 4250 ^5 2c// ^ 6500 and 
—2.0 ^ [Fe/H] ^ 0.5 and testing how the adopted method 
recovers the temperatures and luminosities of the underlying 
synthetic spectra. The accuracy is excellent given that in all 
three infrared bands the IRFM recovers the right tempera- 
ture within 1 K and the bolometric luminosity (i.e. the the- 
oretical value crT^ff) within 0.1% and 0.04% for the ATLAS 
and MARCS models respectively, i.e. at the level of the nu- 
merical accuracy of direct integration of the bolometric flux 
from the spectra. This confirms that the interpolation over 
the grid introduces no systematics and even a poor initial 
estimate of the temperature does not affect the method. 

We expect the dependence of the results on the adopted 
synthetic library to be small, since most of the luminosity 
is actually observed photometrically and the model depen- 
dence for the g-factor and Rtheo (see Section 14.211 is weak 
because we are working in a region of the spectrum largely 
dominated by the continuum. The method is in fact more 
sensitive to the adopted absolute calibration that governs 

Robs • 

The principal contributor to continuous opacity in cool 
stars is due to H~. Blackwell, Lynas-Gray & Petford (1991b) 
have shown that by using more accurate opacities with re- 
spect to previous work, temperatures increased by 1.3% and 
angular diameters decreased up to 2.7%, the effect being 
greatest for cool stars. Even though the dependence of the 
results on the adopted synthetic library is small, the use 
of absolute calibrations derived using the most up-to-date 
model atmospheres (see Appendix A) is of primary impor- 
tance. Though we use a great deal of observational infor- 
mation, the new opacity distribution function (ODF) in the 
adopted grid of model atmospheres is particularly impor- 
tant. With respect to the more accurate but computationally 
time consuming opacity sampling (OS), older ODF models 
can underestimate the IR flux by few percent, translating 
into cooler effective temperatures (Grupp, 2004a). Likewise, 
the better opacities in the UV lead to an increased flux in 
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Figure 9. Comparison between final effective temperatures, bolo- 
metric luminosities and angular diameters for our sample stars, 
when different synthetic libraries (with respect to ATLAS9) are 
used. 



the visual and infrared region. This directly affects the ra- 
tio between the bolometric and monochromatic flux used for 
the Infrared Flux Method (see Section f4.2L so that the lat- 
est model atmosphere have to be preferred (Megessier 1994; 
1998). 

For all our 104 G-K dwarfs we have tested how the 
recovered parameters change by using ATLAS9, MARCS, 
BaSel 2.1 and 3.1 libraries. The comparison in temperatures, 
luminosities and angular diameters is shown in Figure El 

The ATLAS9 and MARCS models show a remarkably 
good agreement through all the temperature range, with dif- 
ferences in the temperature that never exceed a few degrees 
(or 0.2%) and those in the bolometric luminosities and the 
angular diameters being always within 1%. Though not a 
proof of the validity of the models, it is very encouraging 
that the sophisticated physics implemented in independent 
spectral synthesis codes shows such a high degree of consis- 
tency. 

The BaSel models show bigger differences with respect 
to ATLAS9. The BaSel 2.1 model tends to give slightly hot- 
ter temperatures, higher luminosities and smaller angular 
diameters when compared to ATLAS9, though the system- 
atic is oscillating (reflecting the underlying continuum ad- 
justments of the semi-empirical BaSel libraries). 

We conclude by noting that in the range of tempera- 
tures and luminosities studied the dependence of our results 
on the adopted model never exceeds a few percent and for 
the latest models (MARCS and ATLAS9) the agreement is 
well within 1%. Hence we will only present and discuss re- 
sults based on ATLAS9 in the following. 



4.4 Evaluation of the errors: rounding up the 
usual suspects 

Though the IRFM is known to be one of the most accu- 
rate ways of determining stellar temperatures without di- 
rect angular diameter measurements, its dependence on the 
observed photometry and metallicities, on the absolute cali- 
bration adopted and on the library of synthetic spectra make 
the evaluation of the errors not straightforward. In this work 
we have proceeded from first principles and made use of 
high accuracy observations only, thus facilitating the esti- 
mation of errors and possible biases. Since we made use of 
photometric measurements in many bands, we expect the 
dependence on photometric errors to be rather small, small 
random errors in different bands being likely to compensate 
each other. We have also shown in the previous section that 
when the latest model atmospheres are adopted, the model 
dependence is below the one percent level. The effect of the 
absolute calibration can be regarded as a systematic bias; 
once the absolute calibration in different bands is chosen, all 
temperatures scale accordingly. As we show in Appendix A 
the adopted absolute calibration has been thoroughly tested. 



4-4- i Photometric observational error 

The photometric errors in different bands as well as errors 
in metallicity are likely to compensate each other. In or- 
der to check this, for each star we have run 1000 Mon- 
teCarlo simulations, assigning each time random errors in 
BV {RI)c J H K s and [Fe/H]. The errors have been assigned 
with a normal distribution around the observed value and 
a standard deviation of 0.01 mag in BV{RI)c and of 
"h_," "kjnsigcom" in JHKs- Errors in [Fe/H] are those 
given in the papers from which we have collected abun- 
dances. 

The standard deviation of the resulting temperatures 
obtained from the MonteCarlo for each star is typically 
11 K and never exceeds 19 K. Bolometric luminosities have a 
mean relative error of only 0.5% and never exceeding 0.8%; 
for angular diameters the mean error sets to 0.4% and never 
exceeds 0.5%. 



4-4-^ Metallicity and surface gravity dependence 

The interpolation within the grid of synthetic spectra is as a 
function of [M/H]. We have checked the effect of a-elements 
by making the interpolation as a function of [Fe/H] instead, 
finding very little differences in the derived fundamental pa- 
rameters around the solar values. At lower metallicities the 
effect of the a-elements start to be important as expected 
(Figure irnil . The effect is in any case, always within 20 K in 
Tcff and within 2% in bolometric luminosity and 9. 

Spectroscopic metallicities have typical errors well 
within 0.10 dex. As we show in Figure [TTI a change of few 
dex in metallicity has almost no effect on the resulting tem- 
peratures around the solar values, but at lower metallicities 
systematic errors of 0.10-0.20 dex in [Fe/H] can introduce 
biases up to 40 K. This is an important point when using 
large surveys with photometric metallicities that have typi- 
cal errors above 0.10 dex. 

The dependence on the adopted surface gravity is very 
mild for dwarf stars. We have verified that a change of ±0.5 
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Figure 10. Differences between temperatures, bolometric lumi- 
nosities and angular diameters when the [M/H] or [Fe/H] are 
used. 
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Figure 11. Effect of errors in metallicities on the resulting tem- 
peratures. 



dex in \og{g) implies differences that never exceed ~ 30 K in 
temperature and well within 1% for bolometric fluxes and 
angular diameters. 



4-4-3 Systematic error in the absolute calibration 

The errors in the adopted magnitudes and absolute calibra- 
tions of Vega in different bands can be regarded as a sys- 
tematic bias, since once they are selected, the recovered tem- 
peratures, luminosities and angular diameters scale accord- 
ingly. Errors in the observed BV{RI)c magnitudes of Vega 
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Figure 12. Effect of the absolute calibration being systematically 
brighter (circles) or fainter (squares) when uncertainties given in 
Table IX2I correlate. Translated into magnitudes it corresponds to 
a mean shift of Vega's zero-points of ±0.01 mag in the optical 
and ±0.02 mag in the infrared. 



are around 0.01 mag whereas the uncertainties in JHKs 
magnitudes given by Cohen et al. (2003) are within a few 
millimag. As for the random errors on the photometric zero- 
points, it is very likely that uncertainties in Vega's magni- 
tudes compensate each other. 

On the other hand, the uncertainties in the adopted 
absolute calibration mainly come from the uncertainty in 
the ffux of Vega at 5556 A. Since the absolute calibration 
in all bands scales accordingly (even though we have used 
slightly different approaches in optical and infrared, see Ap- 
pendix A), we evaluate the systematic error in temperatures 
for the worst case scenario, i.e. when all errors correlate to 
give systematically higher or lower fluxes. The uncertainties 
adopted are those given in Table El The results are shown 
m Figure m 

Finally, a detailed comparison with data from satellites 
validate within the errors the adopted calibration, though it 
seems to suggest that infrared fluxes should be brightened 
by 1% at most (see Appendix A). If so, the resulting tem- 
peratures would cool down by 10 to 30 K and luminosities 
and angular diameters would increase on average by 0.2% 
and 0.7%, respectively. 



4-4-4 The final error budget 

The primary estimate for errors in T^ff is from the scatter 
in the temperature deduced from J, H and Ks bands that 
reflects photometric errors as well as differences in the corre- 
sponding absolute calibration (see Section [4.211 . As regards 
the uncertainties discussed in Sections 14. 4. 21 and 14.4.11 they 
are of the same order and we quadratically sum them to the 
standard deviation in the resulting J, H and Ks for a fair 
estimate of the global errors on temperatures. 
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Table 2. Mean accuracy of the derived fundamental stellar pa- 
rameters. Notice that either the systematics given here, either 
those shown in Figure lT^ are determined according Table IX2I and 
do not account for possible (and only indicative) shifts in Vega's 
zero-points as a consequence of its rapidly rotating nature (sec 
Appendix A). 



Internal accuracy Systematics 
0.8% 0.3% 
:FboI 0.5% 1.2% 
9 1.7% 0.7% 



5 THE TEMPERATURE SCALE: SOME LIKE 
IT HOT 

In this section we test our IRFM in a number of ways. 
Firstly, we compare its predictions to empirical data for the 
Sun and solar analogs. Secondly, we compare our predicted 
angular diameters to recent measurements with large tele- 
scope interferometers for a small sample of G and K dwarfs. 
Thirdly, we compare our system to other temperature scales 
for G and K dwarfs in the literature. 

Overall, we find that our scale is in good agreement 
with other scales though some puzzles remain. We make 
some suggestions for further work. 



The uncertainties of the bolometric luminosities mostly 
depend on the photometry that accounts for ~ 70-85% of 
the resulting luminosity. The errors from the observations 
have been estimated by summing in quadrature those from 
the MonteCarlo simulation (Section l4.4.H to those due to a 
change of ±0.5 dex in log{g). 

Finally, for the resulting angular diameters we propa- 
gate the errors in temperature and luminosity from equation 
|K| The mean internal accuracy of the resulting temperatures 
is of 0.8%, that of the luminosities of 0.5% and that of the 
angular diameters of 1.7%. 

There are possible systematic errors coming from the 
adopted absolute calibration. As can be seen from Figure IT^ 
such uncertainties give an additional mean error of 0.3% in 
temperature and 1.2% in luminosity that translate into an 
additional systematic error of 0.7% in angular diameters. A 
summary of the uncertainties is given in Table |2| 



4.5 Colour-temperature-metallicity fitting 
formulae 

To reproduce the observed relation Tc// vs. colour and to 
take into account the effects of different chemical compo- 
sitions, the following fitting formula has been adopted (e.g 
Alonso et al. 1996a, Ramirez & Melendez 2005a, Masana et 
al. 2006): 

Ocff = ao + aiX + aaX^ + asXpe/H] 

^-a4[Fe/H]-Ha5[Fe/H]^ (9) 

where 6'c// — 5040/Tcf f , X represents the colour and Ui 
{i = 1,. . . ,5) the coefficients of the fit. In the iterative fitting, 
points departing more than 3cr from the mean fit were dis- 
carded (very few points ever needed to be removed) . All our 
relations were adequately fit by simple polynomials, and we 
did not need to got to higher orders in X to remove possible 
systematics as function of the metallicity (as in Ramirez & 
Melendez, 2005a, who covered a much wider range in tem- 
perature than we do here). 

Figure^]shows the colour-temperature relations in dif- 
ferent bands. The coefficients of the fits, together with the 
number of stars used, the range of applicability and the stan- 
dard deviations are given in Table 01 Note that the very 
small scatter in the relations reflects the high-quality and 
homogeneity of the input data. 



5.1 The Sun and solar analogs 

While the temperature, the luminosity and the radius of the 
Sun are known with great accuracy, its photometric colours 
can only be recovered indirectly. Recently, Holmberg et al. 
(2006) have provided colour estimates for the Sun based 
on those of solar analogs. Another way to obtain colours of 
the Sun levers synthetic photometry, convolving empirical or 
model spectra of the Sun with filter response functions (e.g. 
Bessell et al. 1998) . Both methods lead to estimated colours 
that have typical uncertainties of a few 0.01 mag. Though 
this prevents the use of the Sun as a direct calibrator of 
the temperature scale, it is still useful to compare to its 
estimated colours. 



5.1.1 The colours of the Sun 

We have computed synthetic solar colours from the solar 
reference spectrum of Colina et al. (1996) which combines 
absolute flux measurements (from satellites and from the 
ground) with a model spectrum longward of 9600 A. This 
spectrum, together with those of three solar analogs, is avail- 
able from the CALSPEC library^ . According to Colina et al. 
(1996), the synthetic optical and near- infrared magnitudes 
of the absolutely calibrated solar reference spectrum agree 
with published values to within 0.01-0.03 mag. 

Recently, two new composite solar spectra extending up 
to 24000 A have been assembled by Thuilfier et al. (2004) 
based on the most recent space-based data. The accuracy 
in the UV-visible and near IR is of order 3%. The two so- 
lar spectra correspond to moderately high and moderately 
low solar activity conditions, although the effect of activity 
on the resulting synthetic colours is at or below the mil- 
limag level, except for U ~ B (0.003 mag). Such small differ- 
ences could be due to the solar spectral variability (which 
increases toward shorter wavelengths), the accuracy of the 
measurements or both. In what follows we adopt therefore 
the pragmatic approach of averaging the colours returned 
by the two spectra and we generically refer to the results as 
the Thuilfier et al. (2004) spectrum. 

We have also computed the theoretical magnitudes and 
colours of the Sun predicted by the latest Kurucz and 
MARCS synthetic spectra. Magnitudes and colours com- 
puted from the aforementioned solar spectra are given in 

^ ftp:/ /ftp. stsci.edu/cdbs/cdbs2/calspec/ 
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Figure 13. Empirical colour-temperature relations. Circles are for stars with [Fe/H] ^ —0.5, crosses for stars with [Fe/H] > —0.5. Red 
(cyan) continuous lines are our (Ramirez & Melendez 2005a) fitting formulae at solar metallicity. Dotted lines are for [Fe/H] = —1. In 
the infrared, for stars with T^j f below 5000 K the Ramirez &; Melendez (2005a) calibration shows a larger metallicity dependence than 
indicated by our data. 



Table 3. Coefficients and range of applicability of the colour-temperature relations. 



Colour 


Metallicity range 


Colour range 


ao 


ai 


0.2 


13 


(24 


0.5 


N 


o-(Tc//) 


B - 


V 


[-1.87,0.34] 


[0.419, 1.270] 


0.5121 


0.5934 


-0.0618 


-0.0319 


-0.0294 


-0.0102 


103 


54 


V - 


Rc 


[-1.87,0.34] 


[0.249,0.755] 


0.4313 


1.5150 


-0.7723 


-0.0950 


0.0179 


-0.0033 


101 


55 


{R- 


-I)c 


[-1.87,0.34] 


[0.268, 0.640] 


0.2603 


2.3449 


-1.4897 


-0.1149 


0.0641 


0.0023 


104 


70 


V - 


Ic 


[-1.87,0.34] 


[0.517,1.395] 


0.3711 


0.8994 


-0.2467 


-0.0545 


0.0393 


-0.0010 


101 


49 


V- 


J 


[-1.87,0.34] 


[0.867,2.251] 


0.4613 


0.4118 


-0.0473 


-0.0356 


0.0535 


0.0012 


104 


30 


V - 


H 


[-1.87,0.34] 


[1.144,2.847] 


0.4797 


0.3059 


-0.0252 


-0.0196 


0.0426 


0.0036 


99 


25 


V - 


Ks 


[-1.87,0.34] 


[1.124,3.010] 


0.4609 


0.3069 


-0.0263 


-0.0145 


0.0275 


0.0006 


103 


21 



N is the number of stars employed for the fit after the 3cr clipping and aiT^ff) the final standard deviation of the proposed calibrations. 



Table |1] together with the empirical colours of Holmberg et 
al. (2006). 

Although we have not looked at stellar U magnitudes in 
this study, for completeness with the work of Holmberg et al. 
(2006) we have generated synthetic magnitudes in this band 
too. However, in what follows we focus our discussion from 
B to Ks band (for a discussion of the theoretical solar U — B 
colour see e.g. Grupp 2004b and references therein). We have 
used only Vega to set the zero points. The differences in 
the resulting colours when also Sirius is used for the optical 
bands are given in Appendix A and tend to make the solar 
B — V even redder. Assuming no deficiencies in the synthetic 
spectra, the uncertainties in the derived colours are entirely 
due to the uncertainties in the adopted zero points, i.e. in 
the order of ^ 0.01 mag. 

The CALSPEC and Thuilher et al. (2004) spectra are 
already absolutely calibrated for a distance of 1 AU. For 
the absolute calibration of the Kurucz and MARCS solar 
spectra, we also adopt a distance of 1 AU. The synthetic 
spectra have a temperature Tq = 5777 K and we adopt 
Lq = 3.842 X 10^^ erg s"^ (Bahcall, SereneUi & Basu 2005) 
from which we deduce Rq — ^ Lq /AtkjTq to be used for the 
absolute calibration^. Proceeding this way we can immedi- 

^ The deduced value 9q = 0.00930179 rad compares well with 
that reported in Landolt-Bornstein (1982) Oq = 0.00930484 rad. 



ately compare the recovered temperatures and luminosities 
by means of the IRFM (see Table|SJl with those given above. 
The effect of the absolute calibration of the solar spectra is 
immediately seen in the deduced value of V and reflects in 
the recovered luminosity whereas the temperature depends 
on the colours only. 

An influential direct V measurement of the Sun is 
that of Stebbins & Kron (1957) to which Hayes (1985) 
claimed a further 0.02 mag correction for horizontal extinc- 
tion. Using recent standard V magnitudes for the Stebbins 
& Kron (1957) comparison G dwarfs returns for the Sun 
V = -26.744 ± 0.015, which becomes V = -26.76 ± 0.02 
after accounting for the Hayes correction (Bessell, Castelli 
& Plez 1998). The values deduced using composite and syn- 
thetic spectra are closer to the measurements of Stebbins 
& Kron but fully within the error bars of the Hayes' value. 
Note that through the year the variation of the solar dis- 
tance due to the ellipticity of the Earth's orbit corresponds 
to a V difference of ±0.035 mag. 

In the optical bands synthetic and composite colours 
show a remarkably good agreement with the recent em- 
pirical determination, differences being at most 0.01-0.02 
mag and therefore within the error bars. Synthetic and 
composite spectra also confirm the redder B — V colour 
with respect to the determination of Sekiguchi & Fukugita 
(2000) {B ~V = 0.626) and Rami'rez & Melendez (2005a) 
{B-V = 0.619), though the Thuillier et al. (2004) spectrum 
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Table 4. Magnitudes and colours of the Sun compared to theoretical colours and those deduced using our temperature scale. 



V 


U -B 


B-V 


V-Rc 


(R - l)c 


V - J 


V -H 


V-Ks 


Ref. 




0.173 ±0.064 


0.642 ± 0.016 


0.354 ±0.010 


0.332 ± 0.008 


1.151 ± 0.035 


1.409 ±0.035 


1.505 ±0.041 


(a) 


-26.742 


0.131 


0.648 


0.373 


0.353 


1.165 


1.483 


1.555 


(b) 


-26.743 


0.146 


0.635 


0.374 


0.358 


1.189 


1.528 


1.613 


(c) 


-26.740 


0.101 


0.645 


0.358 


0.349 


1.162 


1.495 


1.556 


(d) 


-26.746 


0.120 


0.658 


0.370 


0.350 


1.166 


1.483 


1.555 


(e) 


-26.753 


0.188 


0.633 


0.360 


0.345 


1.154 


1.485 


1.547 


(f) 






0.651 


0.356 


0.330 


1.150 


1.459 


1.546 


(g) 



(a) Holmberg et al. (2006); (b) Cohna et al. (1996); (c) ThuiUier et al. (2004) (d) ATLAS9 ODFNEW with Grevesse & Sauval solar 
abundances; (e) Kurucz 2004 model with resolving power 100000; (f) MARCS; (g) our temperature scale. 



is slightly bluer than the Neckel & Labs (1984) measure- 
ments used in the Colina et al. (1996) composite spectrum. 
Using the empirical value of B — V = 0.642 in the Ramirez 
& Melendez (2005a) colour-temperature relation implies a 
temperature of the Sun of 5699 K, therefore suggesting the 
need of a hotter temperature scale for this colour. Such an 
offset is of the same order of those found in Section [3.31 for 
optical and infrared bands. 

In the infrared, the only directly measured spectrum is 
that of ThuiUier et al. (2004) . It is much redder than the em- 
pirical determination of Holmberg et al. (2006) and also the 
model spectra. As concerns the other models, J band shows 
good agreement with the empirical determination. Though 
a systematic offset of ~ 0.01 mag appears, it is fully within 
the error bars. The H and Ks bands show larger systematic 
differences in the order of 0.08 and 0.05 mag respectively 
as Figure |S] already suggests. Such large differences are not 
easily understood in terms of model failures. Interestingly, 
an offset of 0.075 mag in the H band zero-point given by 
the 2MASS has also been claimed by Masana et al. (2006) 
to constrain a set of solar analogs to have just the same 
temperature of the Sun. The reason of such a difference re- 
mains unclear as the adopted temperature scale, the model 
and observational uncertainties all might play a role. Among 
other reasons we mention that historically the H filter was 
not defined by Johnson and its values are not easily homog- 
enized, so that further uncertainties might be introduced 
when comparing to various sources, though 2MASS now al- 
lows one to cope with this problem. According to Colina et 
al. (1996) the IR colours deduced from the composite spec- 
trum agree within 0.01 mag with the solar analogs measure- 
ments of Wamsteker (1981) and within 0.03 mag with the 
solar analogs measurements of Campins et al. (1985). The 
synthetic colours obtained using the Colina et al. (1996) 
solar spectrum also agree very closely with those obtained 
when the latest model atmospheres are used. 

Using the colours in Table 2] the temperature of the 
synthetic Kurucz and MARCS spectra are always recovered 
within 7-17 K (see Table differences with the underly- 
ing values likely due to the fact that synthetic solar spec- 
tra are tailored to match the observed abundances, log{g) 
and turbulent velocity whereas the interpolation when ap- 
plying the IRFM is done over a more generic grid. Also the 
temperature of 5802 K deduced when using the synthetic 
magnitudes from Colina et al. (1996) solar spectrum agrees 
very well with the known value of 5777 K whereas the cooler 
temperature returned by the ThuiUier et al. (2004) colours 
is a direct consequence of the much redder IR colours. The 



Table 5. Recovered temperature and luminosity of the Sun. 



Te// (K) 


Luminosity (7^) 


Ref. 


5864 




(a) 


5802 


1.004 


(b) 


5720 


1.008 


(c) 


5770 


0.992 


(d) 


5794 


1.003 


(c) 


5791 


1.005 


(f) 



Table Notes: (a) Holmberg et al. (2006); (b) Colina et al. (1996); 
(c) ThuiUier et al. (2004); (d) ATLAS9 ODFNEW with Grevesse 
& Sauval solar abundances; (e) Kurucz 2004 model with resolving 
power 100000; (f) MARCS. 

temperature of the Sun obtained when using the empirical 
colours of Holmberg et al. (2006) is 87 K hotter. Changing 
the empirical V — H from 1.409 to 1.480 lowers the recov- 
ered temperature to 5811 K; when also V — Ks is reddened 
from 1.505 to 1.550 the recovered temperature then goes to 
5774 K. Therefore - as expected - the temperature recov- 
ered critically depends on the IR magnitudes in H and Ks 
bands. The H band is known to be troublesome in both 
spectral modelling and observations, but a discrepancy in 
the models ~ 0.08 mag is difficult to understand in light of 
the comparison with the observations in Section |3 

The empirical colours of the Sun of Holmberg et al. 
(2006) are deduced interpolating in temperature and metal- 
licity a sample of Sun-like stars with temperatures from 
Ramirez & Melendez (2005b). This scale is some 100 K 
cooler than our own and such a difference can easily ac- 
count for the differences in the IR colours. In comparison, 
the optical colours are almost unaffected by the temperature 
used for the fit since they are primarily dependent on the 
adopted metallicity to fit the solar analogs. 

5.1.2 The solar analogs 

Absolutely calibrated composite spectra of three solar 
analogs, namely P041C, P177D and P330E (Colina & Bohlin 
1997; Bohlin, Dickinson & Calzetti 2001) constitute a set of 
secondary flux standards in addition to the white dwarf stars 
with pure hydrogen atmospheres adopted for the HST UV 
and optical absolute calibration. The three solar analogs also 
determine the absolute flux distribution of the NICMOS fll- 
ters. 

The solar analogs' flux distribution in the ultraviolet 
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and optical regions is based on HST Faint Object Spec- 
trograph (FOS) and Space Telescope Imaging Spectrograph 
(STIS) measurements, whereas longward of 10020 A a scaled 
version of the Colina et al. (1996) absolute flux of the Sun 
is used. 

The STIS solar analogs' flux distribution in the opti- 
cal and longer wavelengths is expected to have uncertainties 
within 2% (Bohlin et al. 2001). Even if the near-infrared 
fluxes of the solar analogs have been constructed from mod- 
els, a thorough comparison with the high-accuracy JHK in- 
frared photometry of Persson et al. (1998) gives confidence 
in the reliability of the adopted spectra. The differences be- 
tween the observed magnitudes and the synthetic ones ob- 
tained by using the spectra of the three solar analogs agree 
within 0.01 mag in J band. The differences in K band are 
less than 0.025 mag whereas in H band the models are from 
0.03 to 0.07 mag brighter. Such differences can partly re- 
flect difficulties in modelling this spectral region; however 
since such a problem is also present when H band synthetic 
magnitudes of the white dwarf calibrators are studied, an al- 
ternative or complementary explanation could simply arise 
from errors in the adopted zero point (Bohlin et al. 2001; 
Dickinson et al. 2002). For the three solar analogs infrared 
magnitudes are also available from 2MASS. The observed 
2MASS and model magnitudes for each of the solar analog 
are given m Table ^ As done by Bohlin et al. (2001) for 
the Persson photometry, the Am values compare the rela- 
tive observed and model magnitudes to the same difference 
for P330E (the primary NICMOS standard): 

Am = [msj„ith(star) — mobs (star)] 

- [ms,,„th(P330E) - mobs(P330E)]. (10) 

The uncertainties in the Am values depend only on the 
repeatability of the photometry and not on uncertainties in 
the absolute calibration of the infrared photometric systems. 

The IR region of the solar analogs' spectra is a scaled 
version of that used for the Sun and the IR synthetic mag- 
nitudes of the solar analogs have been proved to be reliable. 
The colours of the solax analogs are bluer than those ob- 
tained from the solar spectra of Colina, Kurucz and MARCS 
and they are consistent with slightly higher temperatures 
for the solar analogs with respect to the Sun. Nonetheless, 
V - H and V - Ks are still redder by ~ 0.05 and ~ 0.03 
mag respectively and the disagreement is even worse if we 
consider that we are actually comparing solar analogs with 
hotter temperatures (and bluer colours) than the Sun. 

Finally, we test the temperature scale of Section r4.5l via 
11 excellent solar analogs for which accurate B,V, J,H,Ks 
colours are available. Seven are drawn from the top-ten solar 
analogs of Soubiran & Triaud (2004) and 4 from the candi- 
date solar twins of King, Boesgaard & Schuler (2005) . The 
results are shown in Table |5] 

The agreement between our temperatures and those re- 
ported in the two papers is outstanding, the mean difference 
being only ATe// = — 7±50 K our temperatures being only 
slightly cooler. Also, the mean temperature of the 11 solar 
analogs set to 5742 K, thus suggesting that our scale is well 
calibrated at the solar value. 



5.2 Angular diameters 

There are several approaches to deriving effective tempera- 
tures of stars. Except when applied to the Sun, very few of 
them are genuinely direct methods which measure effective 
temperature empirically. 

The direct methods rely on the measurement of the an- 
gular diameter and bolometric flux of the star. These fun- 
damental methods are restricted to a very few dwarf stars, 
although interferometric (e.g. Kervella et al. 2004) and tran- 
sit (Brown et al. 2001) observations have recently increased 
the sample. Nonetheless angular diameters obtained from 
both of the aforementioned methods have to be corrected for 
limb-darkening by using some model. Hence the procedure 
is still partly model-dependent. Observational uncertainties 
stem from systematic effects related to the atmosphere, the 
instrument and the calibrators used (e.g. Mozurkewich et al. 
2003). As to giants and supergiants, where more measure- 
ments are available, the comparison of 41 stars observed on 
both NPOI and Mark III optical interferometers has shown 
an agreement within 0.6%, but with a rms scatter of 4.0% 
(Nordgren, Sudol & Mozurkewich, 2001). 

The different limb-darkened values collected by Kervella 
et al. (2004) for the same stars give an idea of the uncer- 
tainties that might still hinder progress. For dwarf stars in 
particular there are still very few measurements, so that a 
large and statistically meaningful comparison can not yet 
be done. For example, the limb-darkened measurement of 
the dwarf stars HD10700 and HD88230 obtained by Pijpers 
et al. (2003) and Lane et. al (2001), respectively, are ~ 5% 
and ^ 9% smaller than those of Kervella et al. (2004) and 
that are ultimately used by Ramirez & Melendez (2005b) to 
check their scale. Besides, most of the limb-darkening cor- 
rections are done using ID atmospheric models that rely on 
the introduction of adjustable parameters like the mixing 
length. 3D models for the K dwarf a Cen B provide a radius 
smaller by roughly lastat (or 0.3%) compared with what can 
be obtained by ID models (Bigot et al. 2006). However for 
hotter stars the correction due to 3D analysis is expected 
to be larger as a consequence of more efficient convection. 
Interestingly, using parallax to convert the smaller angular 
diameter obtained from 3D models into linear radius returns 
0.863 -R0, which is in better agreement with the asteroseis- 
mic value of 0.857 Rq by Thevenin et al. (2002). For Pro- 
cyon the difference between ID and 3D models amounts to 
roughly 1.6%, implying a correction to Tc// of order 50 K 
(AUende Prieto et al. 2002). 

If temperatures are to be deduced from these diame- 
ters, further dependence on bolometric fluxes, often gath- 
ered from a variety of non-homogeneous determinations in 
literature, enters the game. All these complications eventu- 
ally render 'direct' methods to be not as direct, or model- 
independent, in the end. In any case, direct angular diameter 
measurements are still limited to the nearest (and brightest) 
dwarfs with metallicity around solar. 

In our case, the comparison with direct methods is ham- 
pered since the angular diameter measurements available so 
far are of nearby and bright dwarfs that have unreliable or 
saturated 2MASS photometry. Since we do not have any 
star in our sample with direct angular diameter measure- 
ments, we can only perform an indirect comparison by us- 
ing the stars we have in common with Ramirez & Melendez 
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Table 6. Synthetic and observed 2MASS photometry of HST Solar Analogs. 



J 


J(synth) 


H 


_ff(synth) 


Ks 


Xs(synth) 


Am./ 






Star 


10.864 
12.245 
11.781 


10.869 
12.253 
11.796 


10.592 
11.932 
11.453 


10.552 
11.925 
11.471 


10.526 
11.861 
11.432 


10.479 
11.843 
11.390 


-0.010 
-0.007 


-0.058 
-0.025 


-0.005 
0.024 


P041C 
P177D 
P330E 



P041C = GSC 4413-304; P177D = GSC 3493-432; P330E = GSC 2581-2323. Am values refer to Bohlin's model-2MASS difference 
relative to P330E, as explained in the text. 

Table 7. Magnitudes, colours and temperatures of the Solar Analogs. Temperatures have been deduced using the given colours into the 
IRFM. 



V 


U -B 


B-V 


V-Rc 


(R-I)c 


V - J 


V - H 


V-Ks 


Toff (K) 


Star 


12.005 


0.135 


0.611 


0.346 


0.346 


1.137 


1.454 


1.526 


5828 


P041C 


13.356 


0.139 


0.607 


0.350 


0.347 


1.139 


1.454 


1.528 


5831 


P177D 


12.917 


0.055 


0.604 


0.353 


0.352 


1.148 


1.463 


1.538 


5822 


P330E 



Magnitudes and colours of P177D and P330E have been dereddened with the value of E{B — V) given in Bohlin et al. (2001) and using 
the standard extinction law of O'Donnell (1994) and Cardelli, Clayton & Mathis (1989) in the optical and infrared, respectively. The 
higher temperatures of these solar analogs are in agreement with the bluer colours and the recovered temperature of P041C agree very 
well with the estimated value of 5900 K given in Colina & Bohlin (1997). For the other two solar analogs the reddening is used to adjust 
the continuum in the infrared to match the UV-optical observations thus avoiding any need to estimate the temperature (Bohlin, priv. 
com.) 



(2005b). For the common stars, our angular diameters are 
systematically slightly smaller (~ 3%) than those found by 
Ramirez & Melendez (2005b) which agree well with the di- 
rect determinations of Kervella et al. (2004). Nevertheless, 
for individual measurements all determinations agree within 
the errors. 

Finally, we mention another indirect way to determine 
angular diameters. This is provided by spectrophotomet- 
ric techniques, comparing absolutely calibrated spectra with 
model ones via equation <A3t (Cohen et al. 1996, 1999). 
In the context of the absolutely calibrated spectra assem- 
bled by Cohen and collaborators (see Appendix A) such a 
procedure has shown good agreement with direct angular 
diameter measurements of giants, even though spectropho- 
tometry leads to angular diameters systematically smaller 
by a few percent. Considered the aforementioned scatter for 
angular diameters of giants, the difference is not worrisome. 
However, since the absolute calibration we have used for the 
IRFM is ultimately based on the work of Cohen and collab- 
orators (see Appendix A) , such a difference is interesting in 
light of our results (see also Section f6.H . 



6 COMPARISON WITH OTHER 
TEMPERATURE SCALES 

Most ways to determine effective temperatures are indirect 
methods and to different extent all require the introduc- 
tion of models. Other than via the IRFM, for the spectral 
types covered by this study, effective temperatures may be 
determined via (1) matching observed and synthetic colours 
(Masana et al. 2006), (2) the surface brightness technique 
(e.g. di Benedetto 1998), (3) fitting observed spectra with 
synthetic ones (Valenti & Fischer 2005), (4) fitting of the 
Balmer line profile (e.g. Mashonkina et al. 2003), (5) from 
spectroscopic conditions of excitation equilibrium of Fe lines 



(e.g. Santos et al. 2004, 2005) and (6) line depth ratios (e.g. 
Kovtyukh et al. 2003). 

In what follows we compare our results with those ob- 
tained via these indirect methods. Only single and non- 
variable stars with accurate photometry according to the 
requirements of Section |5| and 12.21 have been used for the 
comparison. 

6.1 Ramirez &: Melendez sample 

The widest application of the IRFM to Pop I and II dwarf 
stars is that of Alonso et al. (1996b). Recently, Ramirez & 
Melendez (2005b) have extended and improved the metallic- 
ities in the sample of Alonso et al. (1996b) recomputing the 
temperatures. As expected, the updated temperatures and 
bolometric luminosities do not significantly differ from the 
original ones, since the absolute infrared fiux calibration (i.e. 
the basic ingredients of the IRFM) as well as the bolometric 
flux calibration used are the same as of Alonso et al. (1994, 
1995). The difference between the old and new temperature 
scale is not significant, though in the new scale dwarf stars 
are some 40 K cooler. 

For 18 stars in common between our study and that of 
Ramirez & Melendez (2005b) we find an average difference 
ATcff — 105 ± 72 K, our scale being hotter. This translates 
into a mean T^ff difference of 2.0%, luminosities brighter 
by 1.4% and angular diameters smaller by 3.3% (Figure [T^ . 
Though not negligible, such differences are within the error 
bars of current temperature determinations. 

The strength of their temperature scale is that the abso- 
lute calibration adopted was derived demanding the IRFM 
angular diameters to be well scaled to the directly measured 
ones for giants (Alonso et al. 1994). The absolute calibra- 
tion of Vega was tuned to return angular diameters that 
matched the observed ones; hence it depends on the input 
angular diameters and IR photometry that were expected to 
be very accurate. Interestingly, the method applied to hot 
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Table 8. Temperature scale applied to Solar Analogs. 







[Fe/H] 


V 


B-V 


J 


H 


Ks 


t:„ (k) 


T.ff (K) 


HD 


168009 


—0.04 


6.30 


0.64(1) 


5.12(0) 


4.83(6) 


4.75(6) 


5801 


5784 


HD 


89269 


—0.23 


6.66 


0.65(3) 


5.39(4) 


5.07(4) 


5.01(2) 


5674 


5619 


HD 


47309 


+0.11 


7.60 


0.67(2) 


6.40(3) 


6.16(1) 


6.08(6) 


5791 


5758 


HD 


42618 


—0.16 


6.85 


0.64(2) 


5.70(1) 


5.38(5) 


5.30(1) 


5714 


5773 


HD 


71148 


-0.02 


6.32 


0.62(4) 


5.15(8) 


4.87(8) 


4.83(0) 


5756 


5822 


HD 


186427 


+0.06 


6.25 


0.66(1) 


4.99(3) 


4.69(5) 


4.65(1) 


5753 


5697 


HD 


10145 


-0.01 


7.70 


0.69(1) 


6.45(0) 


6.11(2) 


6.06(3) 


5673 


5616 


HD 


129357 


-0.02 


7.81 


0.64(2) 


6.58(3) 


6.25(7) 


6.19(2) 


5749 


5686 


HD 


138573 


-0.03 


7.22 


0.65(6) 


6.02(7) 


5.74(2) 


5.66(2) 


5710 


5739 


HD 


142093 


-0.15 


7.31 


0.61(1) 


6.15(8) 


5.91(0) 


5.82(4) 


5841 


5850 


HD 


143436 


0.00 


8.05 


0.64(3) 


6.88(4) 


6.64(9) 


6.54(1) 


5768 


5813 



Solar analogs drawn from Soubiran & Triad (2004) and King, Boesgaard & Schuler (2005). [Fc/H] and T*^ ^ are those reported in the two 
papers. For the stars from Soubiran & Triaud (2004) we have adopted the mean metallicity and effective temperature as given in their 
table 4. T^j f is the final effective temperature obtained by averaging the temperatures given by our (B — V), {V — J), {V — H), {V — Kg) 
calibration. For HD 168009 and HD 186427 {V — J) has not been used because of the poor quality flag associated to J magnitudes. 
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Figure 14. Comparison between effective temperatures, bolo- 
metric luminosities and angular diameters obtained in this work 
(TW) and from Ramirez & Melendez (RM) . We have checked for 
dependencies in all the scales on [Fe/H], Tgff,9 and bolometric 
luminosities and found no significant dependencies. 



(Te/ / > 6000 K) and cool {%f / < 5000 K) stars returned 
absolute calibrations that differed by about 8% in all IR 
photometric bands, thus suggesting that effective temper- 
atures for cool and hot stars and those determined from 
the fRFM were not on the same scale. The cause was some 
sort of systematic error affecting direct angular diameters 
and/or model inaccuracies in predicting fR fluxes; the final 



adopted absolute calibration was then a weighted average of 
that returned for hot and cool stars. 



The comparison of the angular diameters derived by 
Ramirez & Melendez (2005b) for dwarf stars with the 13 
recent interferometric measurement of subgiant and main 
sequence stars (Kervella et al. 2004) seems to imply that 
the adopted absolute calibration holds also in this range 
(Ramirez & Melendez 2004, 2005b). They find good agree- 
ment with direct angular diameter measurements, but one 
should be reminded that a well defined angular diameter 
scale for dwarf stars is not yet available (see Section 15.211 . 
Also, the stars we have in common with Ramirez & Melendez 
(2005b) have all angular diameters 6 < 0.5 mas and differ- 
ences between our and their values lie all below 0.02 mas i.e. 
below the uncertainties of current measurements (see table 
4 of Ramirez & Melendez 2005b). 

We strongly suspect that our hotter temperature scale 
is the result of the 2MASS vs. Alonso absolute calibration as 
the almost constant offset in Figure mi suggests. When we 
apply our fRFM transforming first 2MASS photometry into 
the TCS system used by Alonso, by means of the relations 
given by Ramirez & Melendez (2005b) and adopting the 
TCS filters with the absolute calibration and Vega's zero- 
points in the fR given by Alonso, our temperature scale 
sets onto that of Ramirez & Melendez within 20 K, with 
differences in temperatures, luminosities and angular diam- 
eters well below 1%. The confidence in our adopted absolute 
cahbration and zero-points for Vega comes from the exten- 
sive comparison with ground and space based measurements 
(see Appendix A). Furthermore, we also prefer to avoid any 
transformation from the 2MASS to the TCS system since 
that would introduce further uncertainties of 0.03-0.04 mag 
in the photometry (Ramirez & Melendez, 2005b). 

Finally, we mention another extensive application of the 
fRFM, that of Blackwell & Lynas-Gray (1998). For 10 com- 
mon stars our mean temperatures are hotter by 60 ± 67 K 
(1.0%), our luminosities are brighter by 1.3% and our angu- 
lar diameters smaller by 1.4%. 



© 0000 RAS, MNRAS 000, 000-000 



Accurate fundamental parameters for Lower Main Sequence Stars 21 




4500 



5000 5500 
Teff™ (K) 



6000 



6500 



Figure 15. Same as Figure IT4l but between our work (TW) and 
the Masana, Jordi & Ribas scale (MJR). 
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Figure 16. Same as Figure IT4l but between our work (TW) and 
the Di Benedetto scale (diB). 



6.2 Masana et al. sample 

Recently Masana et al. (2006) have derived stellar effective 
temperatures and bolometric corrections by fitting V and 
2MASS IR photometry. They calibrate their scale by requir- 
ing a set of 50 solar analogs drawn from Cayrel de Strobel 
(1996) to have on average the same temperature as the Sun. 
As a result, they find significant shifts in the 2MASS zero- 
points given by Cohen et al. (2003). Remarkably, the shift 
in H band is found to be 0.075 mag, i.e. of the same order 
of the discrepancy between synthetic and empirical colour 
(once the Ramirez & Melendez scale is adopted) found in 
this band for the Sun (see Section [5. H . 

For 64 common stars our Te/ / scale is cooler by 50±80 K 
(1.0%), our luminosities are fainter by 6.9% and our angular 
diameters smaller by 1.5% (see Figure ITSll . 

Masana et al. (2006) found agreement within 0.3% be- 
tween their angular diameters and the uniform disk mea- 
surements collected in the CHARM2 catalogue of Richichi & 
Percheron (2005), with a standard deviation of 4.6%. Since 
the standard deviation is of the same order of the correction 
between uniform to limb darkened disk, for the comparison 
they adopt the questionable choice of using the uniform disk 
rather than the proper limb-darkened one. 

From 385 stars in common with Ramirez & Melendez 
(2005b), Masana et al. (2006) found that their temperature 
scale is on average 58 K hotter than Ramirez & Melendez. 
From stars in common, we find our temperatures are on av- 
erage 50 K cooler than that of Masana et al. (2006), leading 
us to expect we would be on the same scale as Ramirez & 
Melendez (2005b), yet from stars in common (different stars) 



we find we are 105 K hotter than Ramirez & Melendez, a 
very puzzling result! We worked very hard to impute this in- 
consistency to one or other of the scales, but were unable to 
resolve the problem. This section demonstrates that temper- 
atures, luminosities and angular diameters calibrations for 
our stars, retain systematics of the order of a few percent, 
despite the high internal accuracy of the data. 

6.3 Di Benedetto sample 

Another way to determine effective temperatures and angu- 
lar diameters of stars is via the surface brightness technique. 
We have 9 stars in common with the extensive work of Di 
Benedetto (1998). On average our stars are 50 ± 50 K hot- 
ter, the luminosities are brighter by 0.9% and the angular 
diameters smaller by 1.4% (see Figure [TCI . 

Comparing to the angular diameters from the IRFM 
sample of Blackwell & Lynas Gray (1998), Di Benedetto 
(1998) found in the F-G-K spectral range an overall agree- 
ment well within 1% though with an intrinsic scatter as large 
as 2%. 

6.4 Valenti & Fischer sample 

Valenti & Fischer (2005) have presented a uniform cata- 
logue of stellar properties for 1040 nearby F, G and K stars 
that have been observed by the Keck, Lick and AAT planet 
search programs. Fitting observed echelle spectra with syn- 
thetic ones, they have obtained effective temperatures, sur- 
face gravities and abundances for every star. We have 47 
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stars in common. Except for two of them (HD22879 and 
HD 193901) which depart from the comparison by 242 K and 
497 K respectively, we obtain an excellent average agreement 
ATc// = 6 ± 60 K. We could not single out a reason of such 
a large discrepancy with the two outliers; however we are 
confident in our values since they agree within 80-150 K 
with the temperatures of Ramirez & Melendez (2005b) and 
Masana et al. (2006). 

6.5 Santos et al. sample 

Santos et al. (2004, 2005) have carried out a detailed spec- 
troscopic analysis for 119 planetary-host stars and 95 single 
stars. They obtained spectroscopic temperatures based on 
the analysis of several Fe I and Fe II lines. The comparison 
done by Ramirez & Melendez (2004, 2005b) seems to im- 
ply that such temperatures are hotter by about ~ 100 K. 
However for stars with direct angular diameter measure- 
ments Santos (2005) found an agreement within 7 K with 
the direct measured temperatures reported in Ramirez & 
Melendez (2004). 

After removing a star that departs 280 K (HD 142709), 
the difference between our and Santos Tc// scale is to 
15 ± 81 K, our scale being cooler. The reason of the large 
departure for the one outlier is not clear, but our value agree 
better with Masana et al. (2006). 

Another study with temperatures derived using strictly 
spectroscopic criteria (based on the excitation equilibrium 
of the iron lines) is that of Luck & Heiter (2005). We have 
only 5 stars in common, however the average agreement is 
good (ATc// = 11 ± 142 K). Even though such a result is 
dominated by the scatter, it further suggests that our IRFM 
temperatures agree very well with spectroscopic determina- 
tions. 

6.6 Mashonkina et al. sample 

Balmer line profile fitting allows a very precise deter- 
mination of stellar effective temperature for cool stars 
(Fuhrmann, Axer & Gehren 1993, 1994). Mashonkina et al. 
(2000, 2001, 2003) have extensively used such a technique 
to derive effective temperatures to be used for their detailed 
abundance analysis. We have 9 stars in common and the 
mean difference is 37 ±64 K, our temperatures being hotter. 

6.7 Kovtyukh et al. sample 

The line depth ratios technique (Gray 1989, 1994; Gray & 
Johanson 1991) applied to high S/N echelle spectra is capa- 
ble of achieving an internal precision as high as 5 to 10 K. 
Kovtyukh et al. (2003) has applied such a technique to a set 
of 181 F-K dwarfs and adjusted the zero-point of the scale on 
the solar reflected spectra taken with ELODIE, leading to 
the uncertainty in the zero-point of order 1 K. Unfortunately 
we only have three stars in common (in the range 5000 to 
5500 K), too few to draw any meaningful conclusions. Nev- 
ertheless we find good agreement with a mean difference of 
37 ± 43 K, our stars being cooler. The scale of Kovtyukh 
is in fact hotter by about 90 K when compared with the 
stars in the same range from Ramirez & Melendez (2005b). 
Since the zero-point of Kovtyukh et al. (2003) is calibrated 



on the Sun, one of their conclusion is that either the IRFM 
of Alonso et al. (1996b) and Blackwell & Lynas-Gray (1998) 
or the surface brightness technique of di Benedetto (1998) 
could actually predict a too low temperature for the Sun and 
the solar type stars. As we have seen in Section [5. ll a hotter 
temperature scale for the solar analogs could indeed solve 
the discrepancies between observed and computed colours 
in the infrared. 

6.8 Summary : systematic error remains the 
problem 

Our temperature scale is slightly hotter than other scales 
based on the IRFM, typically by 50 to 100 K. Our scale 
agrees closely to the temperature scale established via spec- 
troscopic methods. The strength of our work is that our 
IRFM was done completely from first principles (from the 
multiband photometry and an adopted absolute flux calibra- 
tion) with the best quality data available. The bolometric 
fluxes of the stars are close to completely observationally 
constrained (only 15 to 30 % of the stellar flux lies outside 
our BV{RI)cJHKs filters). Despite this effort, comparison 
with many other temperature scales forces us to conclude 
that external uncertainties in the temperature scale of lower 
main sequence stars remains of of order ±100 K. The exter- 
nal uncertainty dominates our high internal error, which is 
of order only 40 K. 

A summary of the comparisons of our scale with the 
others mentioned above is given in Table |5] 

7 EMPIRICAL BOLOMETRIC CORRECTIONS 

In this section we derive bolometric corrections for our stars. 
The definition of apparent bolometric magnitude is: 

rriBoi = -2.5 log( J"Boi) ± constant (11) 

where J-boI is the bolometric flux received on the Earth 
as defined in eq. ^ . The usual definition of bolometric cor- 
rection in a given band 

BCc = niBoi -rrn = MboI - (12) 

where BC is to be added to the magnitude in a given 
band to yield the bolometric magnitude. 

Although the definition of bolometric magnitude is 
straightforward, there can be some confusion resulting from 
the choice of zero-point. Originally bolometric corrections 
where defined for the V band only and it had been accepted 
that BCv should be negative for all stars. From the spectral 
energy distribution we expect the bolometric correction in 
the V band to be largest for very hot and for very cool stars. 
The minimum of BCv then occurs around spectral type F, 
which then set the zero-point of the bolometric magnitudes. 
This implied a value BCv for the Sun of between —0.11 mag 
(AUer, 1963) and -0.07 (Morton & Adams, 1968). However 
with the publication of a larger grid of model atmospheres, 
the smallest bolometric correction in Kurucz's grid (1979) 
implied BCv = —0.194. The original choice of the zero-point 
for the bolometric corrections has thus proved troublesome. 
A better method is to adopt a fixed zero point. 

We define the absolute bolometric magnitude of the Sun 
to be MboI,q = 4.74 in accordance with Bessell et al. (1998). 
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Table 9. Comparison with other works. 



Ramirez & Blackwell &; Masana Di Benedetto Valenti & Santos Luck & Mashonkina Kovtyukh 
Melendez Lynas-Gray et ah Fischer ct ah Heiter et ah at ah 



Method 


(a) 


(a) 


(b) 


(c) 


(d) 


(e) 


(e) 


(f) 


(g) 


Stars in common 


18 


10 


64 


9 


45 (2) 


14 (1) 


5 


9 


3 


ATe// (K) 


105 


60 


-50 


50 


6 


-15 


-11 


37 


-37 


a (K) 


72 


67 


80 


50 


60 


81 


142 


64 


43 


AL (%) 


1.4 


1.3 


-6.9% 


0.9% 












Ae (%) 


-3.3 


-1.4 


-1.5 


-1.4 













(a) = IRFM; (b) = Spectral Energy Distribution fit; (c) = Surface Brightness technique; (d) = fitting observed spectra with synthetic 
ones; (e) = excitation equihbrium of Fe lines; (f) = fitting of the Balmer line profile; (g) = line depth ratios. The number of stars in the 
bracket refers to those not counted because of large departures in temperature. Differences are computed our work— others. 



By adopting its measured apparent magnitude Vq = —26.76 
(see Section [5. Ill , the absolute Mv magnitude is thus 4.81 
and the V bolometric correction is then BCv ~ 4.74—4.81 = 
-0.07. 

The absolute bolometric magnitude for a star with lu- 
minosity L, radius R and effective temperature Tc// is 



MboI = -2.5 log (L/Lq) + MboI.q 

= -2.5 log (R'T^ii/RlT^ff^Q) + MboLq 



(13) 



where L = 4nR^aT^j-j-- and L© = 3.842 x 10^^ erg s"^ 

(Bahcall et al. 2005). 

Using eq. p3|l and the definition of absolute magnitude 
= rriQ + 51og(7r) — 10, where vr is the parallax (in mas), 

the bolometric correction M21 is 

BCc = -5 log ^ - 10 log + MboI,q - mc + 10 



= -51og/C 



6/2 



10 log 



c//,0 



+ A/soi,©-mc + 10(14) 



where IC is the conversion factor for the proper unit 
transformation, 6 = 2Rtt/IC is the angular diameter (in mas) 
determined via IRFM and i?© is the solar radius in cm as 
given in Section [5. II From grids of model atmospheres bolo- 
metric corrections can be calculated in a similar manner, 
where the dependence on the radius of the star eliminates 
in the difference between the absolute bolometric and in- 
band magnitudes (e.g. Girardi et al. 2002). 

The comparison between the empirical bolometric cor- 
rections in different bands computed via equation 1141 and 
those predicted by model atmospheres (ATLAS9) is shown 
in Figure [TTl 

This plot can be regarded as the theoretical counter- 
part of Figure |21 Again, the agreement between model at- 
mospheres and empirical data is very good and it should be 
remembered that the model dependence in deducing empir- 
ical bolometric correction from our implementation of the 
IRFM is small (only few 10 %, see SectiongJ. 

A complementary way of deriving stellar integrated 
fluxes is via photometric indices. The integrated flux of a 
star, .T^soi (Earth), depends primarily on its apparent bright- 
ness (especially in Rc and Ic bands, see later), which may 
be measured by its magnitude in different bands. Of lesser 
importance is its temperature, which is a function of a colour 
index and metallicity, ^(X, [Fe/H]). Following Blackwell & 
Petford (1991a), we expect a relation of the form 



The function (?i(X, [Fe/H]) is illustrated in Figure^] in 
which the reduced flux in ^-band, J^soj (Earth) lO'^ "'™? , is 
plotted against different colour indices. For the integrated 
stellar flux we have fltted expressions of the form: 



J^so; (Earth) = W 



bo + biX + biX^ + b-iX"" 



+ fo4X[Fe/H] + fog [Fe/H] + fee [Fe/H]' 



.(16) 



We tried different fltting formulae but this form proved 
the most satisfactory. There are a few flts of this type in 
the literature (Blackwell & Petford 1991a for a subset of the 
colours used here, and Alonso et al. 1995 for the infrared). 
Blackwell et al's flt has no metallicity dependence because 
metallicities were not available. Alonso et al's formula gives 
flts which are good in the IR but do not reproduce the ob- 
served trend in the optical bands. Our fltting formulae work 
well from optical to IR. As for the temperature-colour flt- 
ting, at every iteration points departing more than 3(t from 
the mean flt were discarded. The coefficients of the flts, to- 
gether with the number of stars used, the range of applica- 
bility and the standard deviation of the differences between 
the measured fluxes and the fluxes calculated from the flt- 
ting formula are given in Table [Till 

The scatter in Figure [TSl of the reduced flux in Rc and 
Ic bands is only apparently large, due to the small range 
(in vertical scale) covered by the reduced fluxes in these 
two colours (compare also with Figure I17t . For the sake 
of the proposed calibration, the small range covered by the 
reduced flux in these two bands allows for a very accurate (at 
1% level or below) calibration. Interestingly, for these two 
bands tight relations also exist between bolometric fluxes 
and magnitudes, as seen in Figure IT^ 



8 THE ANGULAR-DIAMETER-COLOUR 
RELATION 

Limb-darkened angular diameters can then be readily ob- 
tained via equation from the temperature and bolomet- 
ric flux calibrations given in Section |1| and [7| In particular, 
when using J magnitudes with various colour indices, we 
have found very tight and simple relations, as can be appre- 
ciated from Figure 1201 Analogous relations in other bands 
are signiflcantly less tight. 

We have fltted relations of the form 



TBoi{^&rt\i) 



lQ-OA,^i^(X, [Fe/H]; 



(15) 



Co 4- ci ^/(i>{mJ^JC) 



(17) 
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Table 10. Coefficients and range of applicability of the absolute flux calibrations. All calibrations are valid for the metallicity range 
-1.87 ^ [Fe/H] ^ 0.34. 

Reduced Colour Colour 6o bi 62 63 64 65 be N cr(%) 

flux range 
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-V 


[0 


419 


1.270] 


1.6521 


8.6205 


-13.7586 


12.5042 


-0.7734 


-0.0644 


-0.0882 


100 


1.7 


V - 


Rc 


[0 


249 


0.755] 


0.4475 


19.0419 


-32.6586 


44.0700 


-0.8436 


0.4652 


0.1417 


98 


2.3 


(R- 


I)c 


[0 


268 


0.640] 


11.3466 


-69.2615 


184.7907 


-105.3337 


4.8308 


-1.1543 


0.1516 


101 


3.5 


V - 


Ic 


[0 


517 


1.395] 


3.9334 


-4.2082 


7.7254 


0.7865 


0.5466 


-0.1401 


0.1045 


99 


1.8 


V - 


- J 


[0 


867 


2.251] 


-0.5405 


9.0544 


-6.6975 


2.5713 


0.5042 


-0.2272 


0.1393 


99 


2.5 
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Z.o4 1 \ 
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-1.0000 
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V - 


Ks 


[1 


124 


3.010] 


-0.3850 


6.4601 


-3.5378 


1.0324 


0.4693 


-0.4049 


0.1425 


102 


2.0 


B - 


- V 


[0 


419 


1.270] 


2.9901 


-1.9408 


1.6146 


0.9227 


-0.0857 


-0.3010 


-0.0800 


96 


1.5 


V- 


Rc 


[0 


249 


0.755] 


2.4783 


-0.6395 


2.4817 


3.3947 


-0.0222 


-0.1769 


-0.0100 


103 


1.2 


(R- 


I)c 


[0 


268 


0.640] 


6.0078 


-28.2448 


69.7676 


-41.6115 


1.2198 


-0.5151 


0.0052 


99 


1.7 


V - 


Ic 


[0 


517 


1.395] 


3.7933 
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N is the number of stars used for the fit after the 3(T clipping and <t(%) the standard deviation of the percentage differences between the 
measured fluxes and the fluxes calculated from the fitting formula. The coefficients of the calibrations bj are given in units of 10~® erg 
cm~2 s~^. 
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Figure 17. Bolomctric corrections calculated from Castclli & Kurucz (2003) ATLAS9 are represented as function ofT^ff. [M/H] equal 
to +0.5 (black line), +0.0 (red line), —0.5 (cyan line), —1.0 (green line), —1.5 (blue line), —2.0 (yellow line). Points correspond to our 
IRFM bolometric corrections for the sample stars in the range [M/H] > 0.25 (black), -0.25 < [M/H] ^ 0.25 (red), -0.75 < [M/H] 
^ -0.25 (cyan), -1.25 < [M/H] ^ -0.75 (green), -1.75 < [M/H] ^ -1.25 (blue). The metallicities given for the model are solar-scaled, 
whereas [M/H] for the stars has been computed using eq. Q- 



where 

(j}{mj,X) = W-° X (18) 

for a given colour index X. 

As can be appreciated from Table 1111 these relations 
show remarkably small scatter (at the few percent level). 
In particular, the upper range covered in angular diameters 
can be used to build a network of small calibrators for fu- 
ture long-baseline interferometric measurements from read- 
ily available broadband photometry. 

The calibration in J — Ks colour is particularly appeal- 
ing, since in this photometric index the effect of interstellar 
extinction is negligible. 

The given angular diameter scale can be tested with the 
G dwarf HD 209458 A for which the linear radius has been 
calculated from planetary transit (Brown et al. 2001). We 
have taken V magnitude from Hipparcos and J, H and Ks 
from 2MASS. Our predicted value of 61 = 0.227±0.003 mas is 
in excellent agreement with that of 9 — 0.228±0.004 mas ob- 
tained by Kervella et al. (2004) using surface brightness re- 
lations calibrated by interferometry. The difference between 



the two values is only 0.4% and therefore well below the 
possible ~ 3% offset previously discussed. Unfortunately, 
both this comparison and those in Section |K| and 16.11 have 
their own limitations since they are indirect, or sensitive 
to photometric uncertainties (for this star, Kervella has to 
convert the 2MASS Ks magnitude into the Johnson system 
used to fit his relations). When our angular size is trans- 
lated into linear radius via Hipparcos parallaxes we obtain 
R — 1.150 ± 0.056 Rq slightly smaller than that of Kervella 
et al. (2004) but in better agreement with the direct estimate 
of i? = 1.146 ± 0.050 i?o obtained by Brown et al. (2001) 
via HST time-series photometry (although all determina- 
tions agree well within the error budget). Unfortunately, in 
this case the uncertainty in parallax dominates the bulk of 
the ~ 5% uncertainty on the radius. 

9 CONCLUSIONS 

We have used the Infrared Flux Method to deduce funda- 
mental stellar parameters of effective temperature, bolomet- 



© 0000 RAS, MNRAS 000, 000-000 



26 Casagrande, Portinari and Flynn 



.5x10-'' 



^ 1.0x10-" — 



u 

■p 5.0x10-5 



> 5.0x10-5 — 



4.0x10-5 — 



"g 2.9x10-5 



2.4x10-5 

X 

2.2x10-5 



^ 2.0x10-5 



^ 1.7x10-5 

X 

1.5x10-5 

o 

K 1.3x10-5 




1 XlO-5 



9.0x10-6 — 



"S 7.0x10-6 



.E 9.0x10-6 

X 

^ 7.0x10-6 
o 

"o 4.9x10-6 
cu 

^ 9.5x10-6 
_c 

X 

6.5x10-6 
I 3.5x10-6 



+ + 

t ++ 



+ + 




13- -KEli 




0t:P -tfl ■ 



+ (R-i)cn- 

V-RcD 

B-vn 

V-l 

v-j □ 

V-H □ 

V-Ks 



_4+ + 



+ + 





rein 



□ □□ 
□ □ □ m 

+ □ □ [gti- □□TTi 
"•" -1- 



■r + 




"""^ ^^^^ ^ 



+ + + 



u □ 




3D-. 



+ + 




■44% +++ 



0.5 



1.0 



1.5 2.0 
Colour index 



2.5 



3.0 



Figure 18. Reduced flux in different f-bands (Earth) 10 '"'^ ) plotted as function of different colour indices. Squares are for stars 

with [Fc/H] < -0.5, crosses for stars with [Fe/H] > -0.5. 
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Figure 19. Left panel: relation between apparent magnitudes in f = Rq (circles) and ( = Ic (crosses) band and bolometric fluxes on 
the Earth (in erg cm~^ s~^). Right panel: relation between absolute magnitudes in the same Rc and Iq bands and absolute bolometric 
magnitudes. The fit is in the form 7 = A + B^, 0"(%) is the standard deviation of the percentage differences between the measured values 
and those calculated from the fit and N is the number of stars employed for the fitting after the 3(t clipping. 



Table 11. Coefficients and range of applicability of the angular diameter calibrations. 



Colour 


Colour range 


Co 


Cl 


TV 


a(%) 


cr(mas) 


B - 


- V 


[0.419, 1.270] 


0.00647 


6.18657 


102 


3.4 


0.008 


V - 


Rc 


[0.249, 0.755] 


0.00340 


8.26244 


103 


2.7 


0.008 


(R- 


I)c 


[0.268, 0.640] 


-0.00385 


8.97758 


103 


2.0 


0.006 


V - 


Ic 


[0.517, 1.395] 


-0.00045 


6.09487 


102 


1.6 


0.005 


V - 


- J 


[0.867, 2.251] 


0.00057 


4.63147 


103 


1.6 


0.005 


V - 


- H 


[1.016, 2.847] 


-0.00022 


4.13607 


102 


1.3 


0.004 


V - 


Ks 


[1.124, 3.010] 


0.00098 


4.01335 


103 


1.2 


0.003 


J - 


Ks 


[0.257, 0.759] 


0.00300 


7.96076 


100 


2.8 


0.007 



A'^ is the number of stars used for the fit after the 3cr clipping. cr(%) is the standard deviation of the percentage differences between the 
angular diameters obtained via IRFM and those calculated from the fitting formula, aijaas) is the standard deviation in terms of mas. 



ric magnitude and angular diameter for a sample of 104 
G and K dwarfs. This semi-direct method is mostly based 
on empirical data obtained from our own observations or 
carefully selected from the literature in order to achieve 
the highest accuracy feasible. Our stars all have excellent 
BV{RI)cJH Ks photometry, excellent parallaxes and good 
metallicities. 

Most of the bolometric flux of our stars is seen in the 
optical and infrared bands covered BV{RI)cJHKs- For 
the remaining 15 to 30% of the flux lying outside these 
bands, model atmospheres have been used. Good to very 
good agreement is found for the colours of the stars and the 
synthetically derived colours of a number of sets of model at- 
mospheres in the literature (ATLAS9, MARCS, BaSel 3.1). 

The zero-points of our temperature, luminosity and an- 
gular diameter scales (Sections |1] and |Hl depend only on 
our adopted absolute calibration of Vega. The advantages 
and disadvantages of our use of Vega are discussed in Ap- 
pendix A. The main advantage is that Vega's absolute cali- 
bration has been extensively studied over a wide wavelength 
range via satellite measurements. This allows us to put firm 



constraints on the systematic uncertainties. A likely disad- 
vantage is that Vega is known to be pole-on and a fast ro- 
tator which prevents the use of a unique Vega model atmo- 
sphere in fitting the observations. 

Our calibration of effective temperature is found to be 
some 100 K hotter than similar determinations in the liter- 
ature, while it is also found to be in close agreement with 
effective temperature scales based on spectroscopic meth- 
ods. To reduce our temperature scale by 100 K, the required 
changes to Vega's 2MASS zero-points and/or absolute cal- 
ibration are significantly larger than the known uncertain- 
ties in these quantities. Neither can the model atmospheres 
be blamed, because we find such good agreement between 
the model atmosphere colours and the empirical data. We 
worked hard to impute inconsistency to one or other of the 
scales, but were unable to resolve the issue beyond any rea- 
sonable doubt. Our angular diameters are smaller by about 
3% with respect to other indirect determinations, but they 
seem more in line with the predictions of asteroseismology 
and 3D model atmospheres. We conclude that temperatures, 
luminosities and angular diameters calibrations for lower 
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Figure 20. Empirical angular diameter-colour relations for our sample of stars. 



Table 12. Predicted angular diameter for HD 209458 A. 

V J H Kg dv-j ev^H dv-Ks 

7.65 6.591 6.366 6.308 0.230 0.225 0.224 0.227 

±0.01 ±0.020 ±0.038 ±0.026 ±0.003 

6^—( are the angular diameters obtained when {V — J), {V — H) and {V — Ks) colour indices are used. Since the accuracy of the 
calibration in these bands is very similar (see Table ITTl . the final angular diameter is the average of the three values weighted by the 
photometric errors. 



main sequence G and K dwarfs retain systematics of the 
order of a few percent. 

The high quality and homogeneity of the data produce 
very tight empirical colour-metallicity-temperature and an- 
gular diameter-colour calibrations. In particular, the rela- 



tion between angular diameters and magnitudes in the J 
band is remarkable and it indicates an high sensitivity of 
this band to angular diameters. Since many lower main se- 
quence stars are known to host planets, our relations can 
be used to accurately determine the physical properties of 
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the parent star thus allowing to effectively break the degen- 
eracy between the properties of the stars and that of the 
planets when extremely high precision data are unavailable 
(e.g. Bakos et al. 2005, Sato et al. 2005). 

Future interferometric measurements of angular diam- 
eters would go a long way to addressing the uncertainties 
we have shown as well as to test our findings. In particular, 
the shape of the visibility function in the second lobe would 
directly probe limb-darkening corrections (e.g. Bigot et al. 
2006) for our stars. Likewise, tying the absolute calibration 
and the observed colours to a solar twin would probably re- 
duce many of these uncertainties. At present good candidate 
solar twins are HD 146233 (Porto de Mello & da Silva) and 
HD 98618 (Melendez et al. 2006). Accurate interferometric 
measurements together with high precision photometry for 
a set of nearby solar-like stars would permit to set the abso- 
lute calibration via solar analogs directly (e.g. Campins et 
al. 1985). Besides, extremely high precision multiple band- 
pass studies of transiting extra-solar planets would provide 
direct angular diameter estimates for our stars, and will fur- 
ther test our temperature scale. 



ACKNOWLEDGMENTS 

We are very grateful to Andrei Berdyugin for instruction in 
the use of the remotely operated telescope at La Palma. We 
thank Steve Willner for very useful discussions on the abso- 
lute calibrations and Johan Holmberg and Jorge Melendez 
for many constructive comments and a careful reading of the 
manuscript, as well as the referee for the same kindness. LC 
and CF acknowledge the hospitality of the Research School 
of Astronomy and Astrophysics at Mount Stromlo where 
part of this work was carried out. We are very grateful to 
the Academy of Finland for considerable financial support 
(grants nr. 206055 and 208792). LC acknowledges the sup- 
port of the Magnus Ehrnrooth Foundation and a CIMO fel- 
lowship. LP further acknowledges the support of a EU Marie 
Curie Intra-European Fellowship under contract MEIF-CT- 
2005-010884. This research has made use of the SIMBAD 
database, operated at CDS, Strasbourg, France. This pub- 
lication makes use of data products from the Two Micron 
All Sky Survey, which is a joint project of the University 
of Massachusetts and the Infrared Processing and Analy- 
sis Center/California Institute of Technology, funded by the 
National Aeronautics and Space Administration and the Na- 
tional Science Foundation. 



REFERENCES 

AUende Prieto C, Asplund M., Lopez R.J.G., Lambert 

D.L., 2002, ApJ, 567, 544 
AUende Prieto C, Barklem P.S., Lambert D.L., Cunha K., 

2004, A&A, 420, 183 
AUer L.H., 1963, The Atmosphere of the Sun and Stars, 

N.Y. Ronald Press 
Alonso A., Arribas S., Marti'nez-Roger C, 1994, A&A, 282, 

684 

Alonso A., Arribas S., Martinez-Roger C, 1995, A&A, 297, 
197 



Alonso A., Arribas S., Marti'nez-Roger C, 1996a, A&A, 
313, 873 

Alonso A., Arribas S., Marti'nez-Roger C, 1996b, A&A, 
117, 227 

Anders E., Grevesse N., 1989, Geochim. Cosmochim. Acta, 
53, 197 

Asplund M., 2003, in lAU Joint Discussion, 15, Elemental 
Abundances in Old Stars and Damped Lyman-a Systems, 
p.8 (astro-ph/0310445) 

Aufdenberg J. P. et al., 2006, preprint ( astro-ph /0603327t 

Aumann H.H. et al., 1984, ApJ, 278, 23 

BahcaU J.N., Serenelh A.M., Basu S., 2005, 
preprint ( astro-ph /0511337 ) 

Bai G.S., Zhao G, Chen Y.Q., Shi J.R., Klochkova V.G, 
Panchuk V.E., Qiu H.M., Zhang H.W., 2004, A&A, 425, 
671 

Bakos G.A., et al., 2006, preprint! astro-ph/0603291 ) 
Beirao P., Santos N.C., Israelian G, Mayor M., 2005, A&A, 
438, 251 

Bessell M.S., 1979, PASP, 91, 589 

Bessell M.S., 1990a, A&AS, 83, 357 

Bessell M.S., 1990b, PASP, 102, 1181 

Bessell M.S., 1995, PASP, 107, 672 

Bessell M.S., 2005, ARA&A, 43, 293 

Bessell M.S., Wood P.R., 1984, PASP, 96, 247 

Bessell M.S., Castelh F., Plez B., 1998, A&A, 333, 231 

Bell R.A., Gustafsson B., 1989, MNRAS, 236, 653 

Bigot L., Kervella P., Thevenin F., Segransan D., 2006, 

A&A, 446, 635 
Blackwell D.E., Shallis M.J., 1977, MNRAS, 180, 177 
Blackwell D.E., Petford A.D., 1991a, A&A, 250, 459 
Blackwell D.E., Lynas-Gray A.E., 1994, A&A, 282, 899 
Blackwell D.E., Lynas-Gray A.E., 1998, A&A, 129, 505 
Blackwell D.E., Shallis M.J., Selby M.J., 1979, MNRAS, 

188, 847 

Blackwell D.E., Petford A.D., Shallis M.J., 1980, A&A, 82, 
249 

Blackwell D.E., Lynas-Gray A.E., Petford A.D., 1991b, 

A&A, 245, 567 
Blackwell D.E., Petford A.D., Arribas S., Haddock D.J., 

Selby M.J., 1990, A&A, 232, 396 
Bohlin R.C., 2000, AJ, 120, 437 
Bohlin R.C., Gilhland R.L., 2004, AJ, 127, 3508 
Bohlin R.C., Dickinson M.E., Calzetti D., 2001, AJ, 122, 

2118 

Brewer M-M., Carney B.W., 2006, AJ, 131, 431 

Brown A.G.A., Arenou F., van Leeuwen F., Lindegren L., 

Luri X., 1997, Hipparcos - Venice 1997 Proc, ESA SP- 

402. ESA Publications Division, Noordwijk, p. 63 
Brown T.M., Charbormeau D., Gilhland R.L., Noyes R.W., 

Burrows A., 2001, ApJ, 552, 699 
Campins H., Rieke G.H., Lebofsky M.J., 1985, AJ, 90, 896 
Cardelh J.A., Clayton G.C., Mathis J.S., 1989, ApJ, 345, 

245 

Castelh F., Kurucz R.L., 2003, in Piskunov N., Weiss 
W.W., Gray D.F., eds, Proc. lAU Symp. 210, Mod- 
elling of Stellar Atmospheres. Astron. Soc. Pac, p.A20 
(astro-ph/0405087) 

Castro S., Rich R.M., Grenon M., Barbuy B., McCarthy 
J.K., 1997, AJ, 114, 376 

Cayrel de Strobel G, 1996, A&ARv, 7, 243 

Ciardi D.R., van Belle G.T., Akeson R.L., Thompson R.R., 



© 0000 RAS, MNRAS 000, 000-000 



30 Casagrande, Portinari and Flynn 



Lada E.A., Howell S.B., 2001, ApJ, 559, 1147 
Clementini G., Gratton R.G., Carretta E., Sneden C, 1999, 

MNRAS, 302, 22 
Code A.D., Bless R.C., Davis J., Brown R.H., 1976, ApJ, 

203, 417 

Cohen M., Walker R.G., Witteborn F.C., 1992a, AJ, 104, 
2030 

Cohen M., Walker R.G., Barlow M.J., Deacon J.R., 1992b, 
AJ, 104, 1650 

Cohen M., Witteborn F.C., Carbon D.F., Davies J.K., 
Wooden D.H., Bregman J.D., 1996, AJ, 112, 2274 

Cohen M., Walker R.G., Carter B., Hammersley P., Kidger 
M., Noguchi K., 1999, AJ, 117, 1864 

Cohen M., Wheaton Wm. A., Megeath S.T., 2003, AJ, 126, 
1090 

Colina L., Bohlin R.C., 1994, AJ, 108, 1931 
Colina L., Bohlin R., 1997, AJ, 113, 1138 
Colina L., Bohlin R.C., CasteUi P., 1996, AJ, 112, 307 
Cousins A.W.J. , 1976, Mem. R. Astron. Soc, 81, 25 
Cousins A.W.J., 1978, Mon. Notes Astron. Soc. S. Afr., 37, 
8 

Cousins A.W.J. , 1980, SAAO Circ, 1, 166 
Crawford, 1975, PASP, 87, 481 
Di Benedetto G.P., 1998, A&A, 339, 858 
Di Benedetto G.P., Rabbia Y., 1987, A&A, 188, 114 
Dickinson M., Sosey M., Rieke M., Bohlin R.C., Calzetti 
D., 2002, in Arribas S., Koekemoer A., Whitmore B., eds.. 
The 2002 HST Cahbration Workshop, STScI Bahimore, 
p.231 

Favata P., Micela G., Sciortino S., 1997, A&A, 323, 809 
Fremaux J., Kupka F., Boisson C, Joly M., Tsymbal V., 

2006, A&A, 449, 109 
Fuhrmann K., Axer M., Gehren T., 1993, A&A, 271, 451 
Fuhrmann K., Axer M., Gehren T., 1994, A&A, 285, 585 
Fulbright J. P., 2000, AJ, 120, 1841 

Girardi L., Bertelli G., Bressan A., Chiosi C, Groenewegen 
M.A.T., Marigo P., Salasnich B., Weiss A., 2002, A&A, 
391, 195 

Golay M., 1974, Introduction to Astronomical Photometry, 

(Dordrecht: Reidel) 
Gratton R.G., Carretta E., Claudi R., Lucatello S., Barbi- 

eri M., 2003, A&A, 404, 187 
Gray D.F., 1989, ApJ, 347, 1021 
Gray D.F., 1994, PASP, 106, 1248 
Gray D.F., Johanson H.L., 1991, PASP, 103, 439 
Gray R.O., Corbally C.J., Garrison R.F., McFadden M.T., 

Robinson P.E., 2003, AJ, 126, 2048 
Grevesse N., Sauval A.J., 1998, Space Sci. Rev., 85, 161 
Grupp F., 2004a, A&A, 426, 309 
Grupp P., 2004b, A&A, 420, 289 

Gulliver A.F., Hill G., Adelman S.J., 1994, ApJ, 429, 81 
Gustafsson B., Bell R.A., Eriksson K., Nordlund A., 1975, 
A&A, 42, 407 

Gustafsson B., Edvardsson B., Eriksson K., Mizuno- 
Wiedner M., J0rgensen U.G., Plez B., 2002, in Hubeny I., 
Mihalas D., Werner K., eds., ASP Conf. Ser. Vol. 288, Stel- 
lar Atmosphere Modeling. Astron. Soc. Pac, San Fran- 
cisco, p. 331 

Hauck B., Mermilliod M., 1998, A&AS, 129, 431 

Hayes D.S., 1985, in Hayes D.S., eds., Proc. lAU Symp. 

Ill, Calibration of fundamental stellar quantities, Reidel: 

Dordrecht, p.225 



Holmberg J., Flynn C, Portinari L., 2006, MNRAS, 367, 
449 

Jimenez R., Flynn C, MacDonald J., Gibson B.K., 2003, 

Sci, 299, 1552 
Johnson H.L., Morgan W.W., 1953, ApJ, 117, 313 
Kervella P., Thevenin F., Di Folco E., Segransan D., 2004, 

A&A, 426, 297 
Kilkenny D., van Wyk P., Roberts G., Marang F., Cooper 

D., 1998, MNRAS, 294, 93 
King J.R,. Boesgaard A.M., Schuler S.C., 2005, AJ, 130, 

2318 

Knude J., 1979, A&A, 71, 344 

Kotoneva E., Flynn C, Jimenez R., 2002, MNRAS, 335, 
1147 

Kovtyukh V.V., Soubiran C, Belik S.I., Gorlova N.I., 2003, 

A&A, 411, 559 
Kupka F., Piskunov N., Ryabchikova T.A., Stempels H.C., 

Weiss W.W., 1999, A&AS, 138, 119 
Kurucz R.L., 1979, ApJS, 40, 1 

Kurucz R.L., ATLAS9 Stellar Atmosphere Programs and 
2 km/s grid. Kurucz CD-ROM No. 13. Smithsonian As- 
trophysical Observatory, Cambridge, MA, 18 

Landolt A.U., 1992, AJ, 104, 340 

Landol-Bornstein, 1982, New Series, Gp VI, Vol 2, As- 
tronomy and Astrophys., Subvolume C, Springer, Berlin- 
Heidelberg-New York 

Lane B.F., Boden A.F., Kulkarni S.R., 2001, ApJ, 551, 81 

Lebreton Y., Perrin M-N., Cayrel R., Baglin A., Fernandes 
J., 1999, A&A, 350, 587 

Lejeune T., Cuisinier P., Buser R., 1997, A&AS, 125, 229 

Leroy J.L., 1993a, A&AS, 101, 551 

Leroy J.L., 1993b, A&A, 274, 203 

Luck R.E., Heiter U., 2005, AJ, 129, 1063 

Malagnini M.L., Morossi C, Rossi L., Kurucz R.L., 1986, 
A&A, 162, 140 

Masana E., Jordi C, Ribas I., 2006, 
preprintHastro-ph/0601049 ) 

Mashonkina L., Gehren T., 2000, A&A, 364, 249 

Mashonkina L., Gehren T., 2001, A&A, 376, 232 

Mashonkina L., Gehren T., Travaglio C, Borkova T., 2003, 
A&A, 397, 275 

Megessier C, 1994, A&A, 289, 202 

Megessier C, 1995, A&A, 296, 771 

Megessier C, 1998, in Bedding T.R., Booth A. J., Davis 
J., eds, Proc. lAU Symp. 189, Fundamental Stellar Prop- 
erties: The interaction between Observation and Theory, 
Kluwer Dordrecht, p. 153 

Melendez J., Dodds-Eden K., Robles J.A., 2006, ApJ, 641, 
133 

Menzies J.W., Cousins A.W.J. , Banfield R.M., Laing J.D., 

1989, SAAO Circ, 13, 1 
Mishenina T.V., Kovtyukh V.V., 2001, A&A, 370, 951 
Mishenina T.V., Kovtyukh V.V., Soubiran C, Travaglio 

C, Busso M., 2002, A&A, 396, 189 
Mishenina T.V., Soubiran C, Kovtyukh V.V., Korotin 

S.A., 2004, A&A, 418, 551 
Morton D.C., Adams T.F., 1968, ApJ, 151, 611 
Mozurkewich D. et al., 2003, AJ, 126, 2502 
Neckel H., Labs D., 1984, Sol. Phys., 90, 205 
Nissen P.E., 1994, Rev. Mex. A&A, 29, 129 
Nordgren T.E., Sudol J.J., Mozurkewich D., 2001, AJ, 122, 

2707 



© 0000 RAS, MNRAS 000, 000-000 



Accurate fundamental parameters for Lower Main Sequence Stars 31 



O'Donnell J.E., 1994, ApJ, 422, 158 
Page! B.E.J. , Tautvaisiene, 1995, MNRAS, 276, 505 
Paulson D.B., Sneden C, Cochran W., 2003, AJ, 125, 3185 
Percival S.M., Salaris M., Kilkenny D., 2003, A&A, 400, 
541 

Persson S.E., Murphy D.C., Krzeminski W., Roth M., 

Rieke M.J., 1998, AJ, 116, 2475 
Peterson D.M. at al., 2004, Proc. SPIE, 5491, 65 
Peterson D.M. at al., 2006, preprint! astro-ph/0603520 ) 
Pijpars P.P., Teixaira T.C., Garcia P.J., Cunha M.S., Mon- 

teiro M.J.P.F.G., Christansen-Dalsgaard J., 2003, A&A, 

406, 15 

Plaz B., 2003, in Munari U., eds., ASP Conf. Ser. VoL 298, 

GAIA Spectroscopy, Science and Technology, p. 189 
Plez B., Brett J.M., Nordlund A., 1992, A&A, 256, 551 
Pompeia L., Barbuy B., Grenon M., 2002, ApJ, 566, 845 
Pompeia L., Barbuy B., Grenon M., 2003, ApJ, 592, 1173 
Porto de Mello G.F., da Silva L., 1997, ApJ, 482, 89 
Povada A., Herrera M.A., Allan C., Cordaro G., Lavalley 

C. , 1994, Rev. Max. A&A, 28, 43 

Price S.D., Paxson C., Engelke C., Murdock T.L., 2004, 

AJ, 128, 889 
Ramirez I., Melendez J., 2004, ApJ, 609, 417 
Ramirez I., Melendez J., 2005a, ApJ, 626, 465 
Ramirez I., Melendez J., 2005b, ApJ, 626, 446 
Reid I.N., van Wyk F., Marang F., Roberts G., Kilkenny 

D. , Mahoney S., 2001, MNRAS, 325, 931 
Richichi A., Percharon I., 2005, A&A, 431, 773 
Ridgway S.T., Joyce R.R., White N.M., Wing R.F., 1980, 

ApJ, 235, 126 

Salaris M., Chieffi A., Straniero O., 1993, ApJ, 414, 580 
Santos N.C., Israehan G., Mayor M., 2004, A&A, 415, 1153 
Santos N.C., Israelian G., Mayor M., Bento J. P., Almeida 

P.C., Sousa S.G., Ecuvillon A., 2005, A&A, 437, 1127 
Sato B., et al., 2005, ApJ, 633, 465 
Saxner M., Hammarback G., 1985, A&A, 151, 372 
Schuster W.J., Nissan P.E., 1989, A&A, 221, 65 
Sekiguchi M., Fukugita M., 2000, AJ, 120, 1072 
Serkowski K., Mathewson D.L., Ford V.L., 1975, ApJ, 196, 

261 

Soubiran C., Triaud A., 2004, A&A, 418, 1089 
Stebbins J., Kron G.E., 1957, ApJ, 126, 266 
Su K.Y.L. et al., 2005, ApJ, 628, 487 

Thevenin F., Provost J., Morel P., Berthomieu G., Bouchy 

F., Carrier F., 2002, A&A, 392, 9 
Thoren P., Fehzing S., 2000, A&A, 363, 692 
Thuilliar G., Floyd L., Woods T.N., Cebula R., Hilsenrath 

E. , Herse M., Labs D., Advances in Space Research, 34, 
256 

Tinbergen J., 1982, A&A, 105, 53 

Tokunaga A.T., Vacca W.D., 2005, PASP, 117, 421 

Tomkin J., Lambert D.L., 1999, ApJ, 523, 234 

Valanti J.A., Fischer D.A., 2005, ApJS, 159, 141 

Vargaly J.L., Farraro R.F., Egret D., Koeppen J., 1998, 

A&A, 340, 543 
Wamstaker W., 1981, A&A, 97, 329 

Wastara P., Lejauna T., Buser R., Cuisinier F., Bruzual G., 

2002, A&A, 381, 524 
Wielen R., Dettbarn C., Jahreii3 H., Lenhardt H., Schwan 

H., 1999, A&A, 346, 675 
Yi S., Demarque P., Kim Y.C., Lee Y.W., Ree C.H., Laje- 

une T., Barnes S., 2001, ApJS, 136, 417 



Yong D., Lambert D.L., Ivans I.L., 2003, ApJ, 599, 1357 



APPENDIX A: BROAD-BAND SYNTHETIC 
PHOTOMETRY AND ABSOLUTE 
CALIBRATION 

The choice of the zero-points is of critical importance to set 
synthetic and real photometry on the same scale. Histori- 
cally, the zero-points of the UBV system were defined "in 
terms of unreddened main-sequence stars of class AO . . . with 
an accuracy sufficient to permit the placement of the zero- 
point to about 0.01 mag" (Johnson & Morgan, 1953). The 
pioneering work of Johnson was continued and extended to 
other bands by Johnson himself over the years; in the south- 
ern hemisphere Cousins (1976, 1978, 1980) and astronomers 
at the South African Astronomical Observatory (SAAO) re- 
fined its accuracy and colours range (Menzies et al. 1989; 
Kilkenny et al. 1998). Similar work was carried out in the 
northern hemisphere by Landolt (1992) at the CTIO. Subtle 
differences exist between different works (see Bessell 2005 for 
a review) and especially between the SAAO and the CTIO 
system (Bessell 1995). Nevertheless, the BV{RI)c system is 
nowadays a well defined one whose main characteristics and 
bandpasses are given in Bessell (1990b). 

The basic equation of synthetic photometry"^ in a given 
band reads as follows (e.g. Girardi et al. 2002): 



-2.5 log 



R^ K' nA)rc(A)dA 
^' /^rc(A)dA 



+ ZPc 



(Al) 



where -F(A) is the flux of the synthetic spectrum (given 
in erg cm~^ s~^ A~^)) ?c(A) the transmission curve of the 
filter comprised in the interval {X,^. , A^^ ) and ZP^ the zero- 
point for the ( band. The ratio between the radius R of 
the star and its distance d from us is known as the dilution 
factor and relates to the angular diameter 9 (corrected for 
limb darkening) via: 



R 
d 



(A2) 



where 9 is in radians. We also note that when all quan- 
tities in the logarithmic term of eq. (lAH are known that 
gives the absolute calibration in C^-band (see Table IX2l for 
a star with spectrum F{X). 

While the zero-points of observational photometry are 
defined on an ensemble of well measured stars, for synthetic 
photometry a common choice is to set the zero-points using 
a reference star for which all the physical quantities in eq. 



^ The formulation given is suited for the energy integration that 
characterize traditional systems based on energy-amplifier devices 
like the Johnson-Cousins. Nevertheless, since the 2MASS trans- 
mission curves are already multiplied by A and renormalized (Co- 
hen ot al. 2003), eq. EH can also be used for the 2MASS photo- 
counting integration. Notice that even if our BV{RI)c measure- 
ments have been done with a photo-counting device (CCD), we 
are placing our photometry on a system of standard stars defined 
with the use of photomultiplier and therefore energy integration 
is the most appropriate. 
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lAlL i.e. the spectrum F{X) (usually synthetic), the dilu- 
tion factor and the observed apparent magnitudes or colour 
indices are known in detail. 

This star is usually Vega, for which we have used the 
magnitudes and colours given in Table EH The optical 
colours are from Bessell (1990a), whereas for the 2MASS sys- 
tem they come from Cohen et al. (2003). The latter were de- 
termined post-facto, comparing the observed 2MASS mag- 
nitudes of 33 stars with the values predicted from the ab- 
solutely calibrated templates built by Cohen et al. (2003) 
within the "Cohen- Walker- Witteborn" (1992a) framework. 

In principle we could adopt a unique spectrum for Vega 
and use it to compute zero-points from B to Ks; in this way 
the dependence on the absolute calibration would cancel out 
when computing colour indices, since all colours scale in the 
same way with the adopted absolute calibration. In practice 
we use a slightly different approach using two models, one for 
the absolute calibration in the optical and one in the infrared 
- each consistent with the corresponding data source. 

The 2MASS magnitudes for Vega are deduced in the 
absolutely calibrated system that is ultimately defined on 
the Kurucz spectra for Vega and Sirius used by Cohen et al. 
(1992b). Therefore, for the JHKs bands we have adopted 
the zero-magnitude fluxes provided by Cohen et al. (2003). 
The zero-magnitude fluxes proposed by Cohen for Vega and 
eight of the primary and secondary stars in the calibration 
network of "Cohen- Walker- Witteborn" has been recently 
confirmed in the 4-24 /im range by the Midcourse Space Ex- 
periment (MSX) flux calibration (Price et al. 2004) to be 
accurate around 1% and thus well within the global error of 
1.46% quoted by Cohen, even though it seems that the fluxes 
of Cohen in the infrared should be brightened by 1%. Also, 
Tokunaga & Vacca (2005) have shown the Vega absolute cal- 
ibration of Cohen et al. (1992b) and the model-independent 
measurements of Megessier (1995) to be identical within the 
uncertainties in the range ~ 1-4 ^.ra. 

For the absolute calibration in the optical bands we 
have adopted a Kurucz (2003) synthetic spectrum for Vega 
with Te// = 9550 K, \og{g) = 3.95, [M/H] = -0.5 and mi- 
croturbulent velocity ^ = 2 km s~^. The resolving power 
used has been 500. This model has recently shown excel- 
lent agreement with the Space Telescope Imaging Spectro- 
graph (STIS) flux distribution over the range 0.17-1.01 
(Bohlin & Gilliland 2004). In particular, in the regions of the 
Balmer and Paschen lines the STIS equivalent widths differ 
from the pioneering work of Hayes (1985) (and used by Col- 
ina et al. 1996 to assemble a composite spectrum of Vega) 
but do agree with the predictions of the Kurucz model. 

Once the synthetic spectrum of Vega in the optical 
bands is chosen, we need to scale it in order to match the 
absolute flux of Vega as measured on the Earth at a certain 
fixed wavelength (A). The monochromatic flux given by the 
model i^modGi(A) is simply related to the same monochro- 
matic flux as measured on the Earth -fEarth(A) by 

J'Earth(A) = (^) ' F,^odcl(A) (A3) 

where the dilution factor (R/d)^ is the ratio between the 
radius of Vega and its distance. In principle the dilution fac- 
tor can be deduced from direct measures of Vega's angular 
diameter. In practice, since such measures are more uncer- 
tain than direct measures of the flux, we proceed the other 



way around. Taking the flux value at 5556 A from Megessier 
(1995) _FEarth(5556 A) = 3.46 x 10"" erg cm'^ s'^ A'^ and 
at the same wavelength from the Kurucz (2003) Vega model 
j;nodei(5556 A) = 5.5015572 x lO'' erg cm'^ s'^ A'^ we ob- 
tain {R/df = 6.2891286 x 10"". This value implies and 
angular diameter of 3.272 mas for Vega, which compares 
very well with the interferometric angular diameter mea- 
surements of 3.24 ± 0.07 mas (Code et al. 1976), 3.28 ± 0.06 
mas (Ciardi et al. 2001) and 3.225 ± 0.032 (Mozurkewich et 
al. 2003). 

We chose to use the Megessier flux in accordance with 
the absolute scale adopted for the STIS (Bohlin & Gilliland 
2004) rather than the value of 3.44 x 10"^ erg cm"^ s"^ A"^ 
±1.45% found in Hayes (1985) and used by Cohen to tie his 
absolute calibration in the infrared. A thorough discussion of 
the differences in the infrared is given by Bohlin & Gilliland 
(2004); here we just mention that this choice together with 
the Kurucz (2003) model in the infrared would give a flux 
density of Vega about 2% lower than in Cohen, thus worsen- 
ing the agreement with MSX. On the other hand, the STIS 
absolute calibration in the optical is expected to be better 
than 2% (Bohlin 2000). In this range, the increase of 0.6% 
relative to Hayes translates in absolute fluxes that are higher 
by the same factor and in differences for the BV{RI)c zero- 
points of 0.006 mag. 

Recent work by Gulliver et al. (1994), Peterson et al. 
(2004, 2006) and Aufdenberg et al. (2006) indicate that 
Vega is pole-on and fast rotator. Thus standard model at- 
mospheres are not appropriate for Vega, as discussed by 
Bohlin & Gilliland (2004). Nonetheless STIS in the opti- 
cal and MSX in the infrared have confirmed the adopted 
absolute calibrations to be accurate at the percent level. It 
is interesting to note that Peterson et al. (2006) give an es- 
timate of the shift in Vega's zero-points as a consequence of 
the higher infrared flux due to its rapidly rotating nature. 
Though only indicative, as the authors say, the shifts for 
their latest non-rotating synthetic spectrum in the infrared 
amount to 0.05 — 0.07 mag. The net effect of such shifts 
would be to reduce the temperatures recovered via IRFM. 
Considering that in the infrared we are adopting the abso- 
lute calibration given by Cohen (already some ~ 2% brighter 
as compared to the latest Kurucz model) this reduces the 
size of the required shifts to about 0.02 — 0.04 mag. The 
agreement within 1% with MSX would however require a 
shift of only 0.01 mag. Since the proposed shifts are only 
indicative and the measurements confirm the adopted cali- 
bration within the errors, we adopt the pragmatic choice of 
using the model that at any given range better agrees with 
the data. It is evident from the above discussion that there is 
no consistent published atmospheric model for Vega at both 
visible and infrared wavelengths and this justifies our adop- 
tion of the Kurucz (2003) or Cohen models at different wave- 
lengths. Likewise, the use of composite model atmospheres 
for Vega is also advocated from Peterson et al. (2006). For 
synthetic photometry the use of composite spectra is often 
avoided on the reason that fixing the zero-points, possible 
physical deficiencies in the synthetic spectrum of Vega are 
likely to be present in the synthetic library and therefore 
they would cancel out (e.g. Girardi et al. 2002) . In our case 
the situation is not so clear, since synthetic spectra give now 
excellent comparison with the real ones and the deficiencies 
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Table Al. Observed magnitudes for Vega. 



B 


V 


Rc 


Ic 


J 


H 


Ks 


Ref. 


0.02 


0.03 


0.039 


0.035 


-0.001 


+0.019 


-0.017 


Bessell 1990a 
Cohen et al. 2003 



Table A2. Absolute calibration and effective wavelength of the ground-based optical-IR photometry of Vega used in this work. Quantities 
tabulated correspond to the definition of the zero magnitude in each filter. 



Band 


^eff 


Monochromatic Absolute Flux 


Uncertainty 


Ref. 




A 


erg cm~-^ s^^A~^ 


erg cm~^ s^^A~^ 




B 


4362 


6.310E-9 


6.310E-11 


This paper 


V 


5446 


3.607E-9 


3.607E-11 


This paper 


Rc 


6413 


2.153E-9 


2.153E-11 


This paper 


Ic 


7978 


1.119E-9 


1.119E-11 


This paper 


J 


12285 


3.129E-10 


5.464E-12 


Cohen et al. (2003) 


H 


16385 


1.133E-10 


2.212E-12 


Cohen et al. (2003) 




21521 


4.283E-11 


8.053E-13 


Cohen et al. (2003) 



The effective wavelengths associated to each filter axe computed in accordance with Appendix B. The error in BV{RI)c bands are 
computed assuming an arbitrary uncertainty of 1% to the absolute flux. Lacking the actual uncertainties in the measurement of the 
filters' transmission curve, the adopted uncertainty is consistent with a wavelength— independent fllter uncertainty in addition to the 0.7% 
absolute uncertainty at 5556 A from Megessier (1995). 

for Vega are likely to be due to the particular nature of this 
star. 

Sirius is often used as an additional fundamental colour 
standard (e.g. BcsscU, CastcUi & Plcz 1998), allowing one 
to control the possible problems of putative variability and 
an IR excess for Vega. First seen by Aumann (1984) be- 
yond 12 fim, Vega's variability at shorter wavelengths has 
long been a matter of discussion (see e.g. Bessell, Castelli 
& Plcz 1998; Ciardi et al. 2001), though the rapidly ro- 
tating model could actually resolve the controversy (Peter- 
son at al. 2006). The most recent data from MSX (Price 
et al. 2004) and SPITZER (Su et al. 2005) support mod- 
els where Vega's infrared excess and variability is due to a 
cold (lust disk. The IR excess increase steadily longward of 
12 yum whereas at 4 /xm the differences between model atmo- 
sphere and measurements axe entirely consistent within the 
uncertainties and therefore should not affect the optical and 
near-infrared region (0.36 - 2.4 /xm) we are working with. 
For the sake of completeness we have tested the differences 
in the derived BV{RI)c zero-points when including Sirius in 
the calibration. Adopting the latest Kurucz model available 
for it and the colours from Bessell (1990a), the differences 
in the derived zero-points between the use of Vega only or 
Vega plus Sirius range from 0.008 mag in _B — to 0.002 mag 
in {R — I)c- Such differences axe below the uncertainties in 
magnitude and colours. In addition, since absolute calibra- 
tion and magnitudes for Sirius in the 2MASS system are not 
available, we have decided to adopt only Vega as our stan- 
dard star. This also makes clearer the comparison between 
synthetic colours and the IRFM that strongly depends on 
the Vega absolute calibration. 

The transmission curves used come from Bessell (1990b) 
for the Johnson- Cousins BV(RI)c system and from Cohen 
et al. (2003) for the 2MASS JHKs system. 

The 2MASS transmission curves carefully incorporate 
the effect of the optical system, the detector quantum effi- 



ciency and the site-specific atmospheric transmissions. The 
effect of the telescope optics on the estimated BV{RI)c 
passband are not counted, but in the V band this has been 
proved to cause uncertainty below few millimag (Colina & 
Bohlin, 1994), and thus smaller than any possible observa- 
tional error we are going to compare with. 

We note that all our photometry is reduced to zero air- 
mass and when we refer to measurements on the Earth, we 
always mean at the top of the Earth's atmosphere. 



APPENDIX B: CHARACTERISTIC 
PARAMETERS OF A PHOTOMETRIC SYSTEM 

In the present work we have reduced broadband (heterochro- 
matic) measurement in a band f to a monochromatic flux 
(.Ff) by means of the following relation: 

J^'f F(A)Tc(A)dA 

where the integration is done in the interval A,;^ , Af j, 
that comprise a given passband (A) . We remark that the 
above formulation is suited for both energy integration in 
the optical and photo-counting in the IR, provided that we 
are dealing with the Bessell (1990b) and Cohen et al. (2003) 
transmission curves (see Appendix A). 

Assuming that the function F{X) is continuous and that 
T({X) does not change sign in the interval A,;; , A^^: , the gen- 
eralization of the mean value theorem states that there is at 
least one value of A in the interval A^^ , A^^, such that 

F{Xi) / Tc(A)dA= / F(A)Tf(A)dA. (B2) 

Rearranging this, we obtain 



© 0000 RAS, MNRAS 000, 000-000 



34 Casagrande, Portinari and Flynn 



f, ' F(A)rc(A)dA 

= i^(AO = , (B3) 

/,;/Tc(A)dA 

where A^ is the isophotal wavelength. The isophotal 
wavelength is thus the wavelength which must be given to 
the monochromatic quantity F{Xi) obtained from a hete- 
rochromatic measurement. 

Stellar spectra do not necessarily satisfy the require- 
ments of the mean value theorem for integration, as they 
exhibit discontinuities. Although the mean value of the in- 
trinsic flux is well-defined (and we have extensively used 
it), the determination of the isophotal wavelength becomes 
problematic because spectra contain absorption lines and 
hence the definition can yield multiple solutions. 

Several authors avoid using the isophotal wavelength 
and introduce the effective wavelength defined by the fol- 
lowing expression 

f'' XF{\)T^{\)A\ 

^cjf^^ . (B4) 

1,1^ F(A)rc(A)dA 

This is the mean wavelength of the passband as 
weighted by the energy distribution of the source over the 
band. The effective wavelength is thus an approximation for 
Xi, however since the monochromatic magnitude at the effec- 
tive wavelength is almost identical to that at the isophotal 
wavelength (Golay 1974), we have used the effective wave- 
length through all our work. 

We have verified that the detailed choice of the wave- 
length associated to a monochromatic flux is not crucial for 
the derived correction of the C factor (eq. EJ. 

When dealing with photometric data, a way to obtain 
monochromatic flux from heterochromatic measurements is 
provided by the q factor. Ideally the g-factor (eq. ISJ should 
be determined from spectroscopic data but in practice we 
rely on a grid of models. Since the absolute monochromatic 
flux for a given band is determined as described in Appendix 
A, the definition of the g-factor for a star with spectra F{X) 
is now: 

f^^ T:H(A)dA 

g(AiR) = (B5) 

^/;/f(A)r,.dA 

This definition slightly differs from that given by Alonso 
et al. (1994, 1996b) to refiect the different absolute cali- 
bration we have adopted (i.e. argument of the logarithm 
in equation IA1|I . However, the different definition of the q 
factor has only minor effects on the recovered stellar param- 
eters. Using the definition of Alonso et al. (1996b) would 
go in the direction of making temperatures hotter by 20 
to 30 K on average, the difference increasing with increasing 
temperature. On average also luminosities would be brighter 
by 0.14% and angular diameters smaller by 0.9%. 
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Table 1. Observable and physical quantities for our sample stars. 



name 


HIP 




7r±5 
(mas 


TV 




(K) 




e±Ae 

(mas) 


rriBol 




V 


B 


-V 


V-Rc 


{R-I)c 






J 








H 








[Fe/H] 


[a/Fe] 




HD 


3765 


3206 


57. 


,90 







.98 


5002 


± 


57 


0, 


,426 




0.010 


7.063 


7, 


.360 


0, 


,940 


0.540 


0.440 


5, 


,69 


± 


0.02 


5. 


,27 


± 


0.05 


5. 


,16 


it 


0.02 


0.01 


0.05 


a 


HD 


4256 


3535 


45, 


,43 


± 


0, 


,95 


4966 




46 


0, 


,326 


± 


0.007 


7.675 


7, 


,990 


0, 


,990 


0.560 


0.460 


6, 


,31 


± 


0.02 


5, 


.87 


± 


0.04 


5, 


.74 


it 


0.03 


0.34 


-0.08 


b 


HD 


4308 


3497 


45, 


,76 


± 


0, 


.56 


5781 


± 


88 


0, 


,423 


± 


0.013 


6.449 


6, 


,560 


0, 


,650 


0.373 


0.352 


5, 


,37 


± 


0.02 


5, 


.10 


± 


0.02 


4, 


.95 


± 


0.02 


-0.40 


0.05 


c 


HD 


4747 


3850 


53, 


,09 


± 


1 


.02 


5406 


± 


53 


0, 


,375 


± 


0.008 


7.002 


7. 


,155 


0, 


,790 


0.420 


0.400 


5, 


,81 


± 


0.02 


5, 


.43 


it 


0.05 


5, 


.30 


± 


0.03 


-0.21 


0.06 


d 


HD 


5133 


4148 


71, 


,01 


± 


0, 


,78 


4956 


± 


52 


0, 


,474 


± 


0.011 


6.870 


7, 


,180 


0, 


,940 


0.540 


0.460 


5, 


,54 


± 


0.03 


5, 


.05 


it 


0.03 


4, 


.89 


± 


0.03 


-0.10 


0.04 


e 


HD 


8638 


6607 


24, 


,75 


± 


1 


.22 


5568 


± 


59 


0, 


,207 


± 


0.005 


8.164 


8, 


,300 


0, 


,690 


0.385 


0.365 


7, 


,00 


± 


0.02 


6, 


.68 


it 


0.03 


6, 


.62 


it 


0.02 


-0.50 


0.12 


c 


HD 


10002 


7539 


29. 


,72 


± 


1 


.08 


5258 


± 


63 


0, 


.254 


it 


0.006 


7.968 


8, 


.140 


0, 


.860 


0.450 


0.390 


6, 


,68 


± 


0.02 


6, 


.34 


it 


0.03 


6, 


.21 


it 


0.02 


0.20 


0.00 


b 


HD 


10853 


8275 


43. 


.37 


± 


1, 


.38 


4638 


± 


78 


0, 


.267 


± 


0.009 


8.405 


8, 


.910 


1, 


.040 


0.625 


0.537 


6, 


.93 


± 


0.02 


6, 


.44 


it 


0.03 


6, 


.32 


± 


0.02 


-0.74 


0.21 


f 


HD 


11020 


8346 


25. 


,79 


± 


1, 


.28 


5288 


± 


31 


0, 


.172 


± 


0.002 


8.790 


8, 


.980 


0, 


.810 


0.443 


0.405 


7, 


.55 


± 


0.03 


7, 


.13 


it 


0.05 


7, 


.06 


it 


0.02 


-0.28 


0.21 


b 


HD 


11130 


8543 


37. 


,33 







.97 


5219 




44 


0, 


.270 




0.005 


7.868 


8, 


.070 


0, 


.758 


0.430 


0.403 


6, 


,60 




0.02 


6, 


.20 


it 


0.02 


6, 


.11 


it 


0.02 


-0.57 


0.15 


f 


HD 


11373 


8867 


45, 


16 


± 


1 


,09 


4738 




58 


0, 


.301 


± 


0.008 


8.052 


8. 


,483 


1, 


,020 


0.607 


0.489 


6, 


.61 


± 


0.02 


6, 


,11 


± 


0.04 


6, 


.02 


± 


0.02 


-0.45 


0.12 


f 


HD 


16623 


12364 


15, 


18 




1 


,15 


5858 




55 


0, 


,149 


± 


0.003 


8.657 


8, 


,760 


0, 


,600 


0.335 


0.335 


7, 


.62 


± 


0.02 


7, 


.30 


± 


0.04 


7, 


.28 


± 


0.04 


-0.60 


0.13 


c 


HD 


17190 


12926 


38, 


,95 


± 


1 


,13 


5139 


± 


74 


0, 


,302 


± 


0.009 


7.692 


7, 


,878 


0, 


,827 


0.418 


0.404 


6, 


.37 


d= 


0.02 


6, 


.00 


it 


0.04 


5, 


.87 


it 


0.03 


-0.11 


-0.00 


b 


HD 


18632 


13976 


42, 


,66 


± 


1 


,22 


4984 


± 


87 


0, 


,316 


± 


0.011 


7.726 


8, 


.020 


0, 


,918 


0.529 


0.439 


6, 


.32 


± 


0.02 


5, 


.95 


it 


0.04 


5, 


.84 


it 


0.02 


0.18 


-0.01 


g 


HD 


21197 


15919 


66, 


,15 


± 


0, 


.95 


4577 


± 


60 


0, 


,453 




0.013 


7.314 


7, 


,868 


1, 


,162 


0.684 


0.565 


5, 


.83 


± 


0.02 


5, 


.31 


it 


0.03 


5, 


.12 


it 


0.02 


0.27 


0.13 


j 


HD 


22879 


17147 


41. 


,07 


± 


0, 


.86 


5930 




67 


0, 


.379 




0.009 


6.577 


6, 


.679 


0, 


.550 


0.325 


0.337 


5, 


.59 


± 


0.02 


5, 


.30 


it 


0.03 


5, 


.18 


it 


0.02 


-0.91 


0.34 


b 


HD 


25329 


18915 


54. 


.14 


± 


1, 


.08 


4796 


± 


62 


0, 


.279 


± 


0.008 


8.164 


8, 


.519 


0, 


.876 


0.495 


0.478 


6, 


.77 


± 


0.03 


6, 


.30 


it 


0.04 


6, 


.20 


it 


0.01 


-1.80 


0.49 


i 


HD 


25665 


19422 


54. 


,17 


± 


0, 


.79 


4902 


± 


65 


0, 


.380 


± 


0.011 


7.398 


7, 


.757 


0, 


.873 


0.567 


0.443 


5, 


.99 


± 


0.03 


5, 


.54 


it 


0.03 


5, 


.46 


it 


0.02 


-0.03 


-0.00 


b 


HD 


25673 


19007 


24. 


,23 


± 


1, 


.53 


5132 ± 


113 


0, 


.145 


± 


0.007 


9.291 


9 


.532 


0, 


.823 


0.467 


0.443 


7, 


.96 


± 


0.02 


7, 


.65 


it 


0.04 


7, 


.46 


it 


0.03 


-0.53 


0.15 


i 


HD 


28946 


21272 


37. 


.33 


± 


2, 


.58 


5308 ± 


55 


0, 


.271 


± 


0.006 


7.787 


7, 


.940 


0, 


.764 


0.392 


0.400 


6, 


.52 


± 


0.02 


6, 


.16 


it 


0.03 


6, 


.07 


it 


0.02 


-0.03 


-0.04 


h 


HD 


30501 


22122 


48. 


,90 







.64 


5153 ± 


40 


0, 


.350 




0.006 


7.360 


7, 


.580 


0, 


.880 


0.480 


0.425 


6, 


.06 


it 


0.02 


5, 


.64 


it 


0.02 


5, 


.53 


it 


0.02 


0.13 


-0.02 


j 


HD 


31128 


22632 


15, 


,55 


± 


1 


.20 


6089 ± 


81 


0, 


,118 


± 


0.003 


8.996 


9, 


,135 


0, 


,489 


0.309 


0.325 


8, 


.03 


± 


0.02 


7, 


.80 


it 


0.03 


7, 


.74 


it 


0.02 


-1.54 


0.41 


i 


HD 


34328 


24316 


14, 


,55 


± 


1 


.01 


6025 ± 


70 


0, 


.106 


± 


0.003 


9.275 


9, 


,427 


0, 


,496 


0.312 


0.335 


8, 


.32 


± 


0.02 


8, 


.05 


it 


0.03 


8, 


.00 


it 


0.03 


-1.69 


0.38 


i 


HD 


37792 


26688 


20, 


,61 


± 


0, 


.89 


6556 ± 


100 


0, 


.190 


± 


0.006 


7.641 


7, 


,697 


0, 


,419 


0.249 


0.268 


6, 


.83 


± 


0.02 


6, 


.68 


it 


0.03 


6, 


.57 


it 


0.02 


-0.60 


0.23 


k 


HD 


39715 


27918 


37, 


,57 


± 


1 


.26 


4759 


± 


37 


0, 


.253 


± 


0.004 


8.410 


8, 


.830 


1, 


.010 


0.605 


0.495 


6, 


.99 


it 


0.03 


6, 


.47 


it 


0.04 


6, 


.35 


it 


0.02 


-0.04 


0.01 


b 


HD 


51219 


33382 


31, 


,15 


± 


0, 


.85 


5584 


± 


58 


0, 


.306 


± 


0.007 


7.303 


7, 


.410 


0, 


.699 


0.377 


0.357 


6, 


.15 


it 


0.02 


5, 


.83 


it 


0.03 


5, 


.72 


it 


0.02 


-0.09 


0.04 


1 


HD 


53545 


34146 


15, 


,74 


± 


0, 


.91 


6393 


± 


55 


0, 


.169 


± 


0.003 


8.005 


8, 


.047 


0, 


.456 


0.268 


0.271 


7, 


.16 


it 


0.03 


6, 


.90 


it 


0.03 


6, 


.86 


it 


0.02 


-0.29 


0.11 


i 


HD 


53927 


34414 


44. 


,92 


± 


1, 


.43 


4882 


± 


35 


0, 


.288 


± 


0.005 


8.018 


8, 


.351 


0, 


.889 


0.529 


0.450 


6 


.64 


it 


0.02 


6, 


.14 


it 


0.02 


6, 


.06 


it 


0.03 


-0.37 


0.09 


h 


HD 


57901 


35872 


39. 


.91 


± 


1, 


.14 


4865 


± 


48 


0, 


.311 


± 


0.007 


7.867 


8, 


.224 


0, 


.962 


0.564 


0.452 


6, 


.46 


± 


0.03 


5, 


.98 


± 


0.03 


5, 


.89 


± 


0.02 


0.00 


0.03 


f 


HD 


59374 


36491 


20. 


,00 




1, 


.66 


5839 


± 


58 


0, 


.171 


± 


0.004 


8.373 


8, 


.485 


0, 


.541 


0.324 


0.326 


7, 


.33 


it 


0.02 


7, 


.03 


it 


0.02 


6, 


.98 


it 


0.03 


-0.88 


0.29 


i 


HD 


59468 


36210 


44. 


,43 


± 


0, 


.53 


5625 


± 


57 


0, 


,111 


± 


0.009 


6.630 


6, 


,720 


0, 


,700 


0.377 


0.355 


5, 


.51 


it 


0.03 


5, 


.15 


it 


0.02 


5, 


.04 


it 


0.02 


0.02 


0.00 


b 


HD 


59747 


36704 


50. 


,80 


± 


1, 


.29 


5063 


± 


54 


0, 


,319 


± 


0.008 


7.443 


7, 


,695 


0, 


,877 


0.493 


0.416 


6, 


.09 


± 


0.02 


5, 


.66 


it 


0.02 


5, 


.59 


it 


0.02 


0.06 


-0.04 


b 


HD 


65486 


38939 


54. 


,91 


± 


1, 


.15 


4678 


± 


29 


0, 


,326 


± 


0.005 


7.934 


8, 


,400 


1, 


,050 


0.605 


0.530 


6, 


.51 


± 


0.02 


5, 


.94 


it 


0.02 


5, 


.83 


it 


0.02 


-0.24 


0.07 


e 


HD 


67199 


39342 


57. 


,88 


± 





.58 


5148 


± 


44 


0, 


,425 


± 


0.008 


6.943 


7, 


,166 


0, 


,889 


0.480 


0.416 


5, 


.65 


it 


0.02 


5, 


.23 


it 


0.03 


5, 


.11 


it 


0.02 


0.03 


-0.04 


b 


HD 


73524 


42291 


36. 


,40 


± 





.57 


5978 


± 


60 


0, 


,389 


± 


0.009 


6.486 


6, 


.530 


0, 


,600 


0.330 


0.310 


5, 


.46 


it 


0.02 


5, 


.21 


it 


0.03 


5, 


.10 


it 


0.03 


0.12 


-0.00 


b 


HD 


85512 


48331 


89. 


,67 


± 


0, 


.82 


4429 ± 


134 


0, 


.557 


± 


0.034 


7.008 


7, 


.670 


1, 


.190 


0.720 


0.622 


5, 


.45 


it 


0.02 


5, 


.00 


it 


0.03 


4, 


.72 


it 


0.02 


-0.18 


0.06 


d 


HD 


87007 


49161 


23. 


,23 


± 


1, 


.41 


5260 ± 


44 


0, 


.187 


± 


0.004 


8.632 


8, 


.802 


0, 


.821 


0.445 


0.391 


7, 


.35 


it 


0.02 


6, 


.97 


it 


0.03 


6, 


.89 


it 


0.04 


0.27 


0.00 


j 


HD 


89813 


50782 


37. 


.30 


± 


1, 


.31 


5372 ± 


51 


0, 


.286 


± 


0.006 


7.618 


7, 


.769 


0, 


.784 


0.423 


0.383 


6, 


.40 


± 


0.03 


6, 


.03 


± 


0.04 


5, 


.91 


± 


0.02 


-0.17 


0.06 


f 


HD 


90663 


51254 


39. 


,75 


± 


1, 


.54 


4936 ± 


43 


0, 


.258 


± 


0.005 


8.209 


8, 


.530 


0, 


.920 


0.541 


0.439 


6, 


.85 


it 


0.02 


6, 


.39 


it 


0.02 


6, 


.26 


it 


0.02 


-0.29 


0.08 


f 


HD 


92786 


52470 


39, 


,56 


± 





.99 


5301 ± 


69 


0, 


.265 


± 


0.007 


7.841 


8, 


,020 


0, 


,750 


0.427 


0.385 


6, 


.57 


± 


0.02 


6, 


.23 


± 


0.02 


6, 


.13 


it 


0.02 


-0.29 


0.08 


f 


HD 


94028 


53070 


19, 


,23 


± 


1, 


.13 


6091 ± 


68 


0, 


.180 


± 


0.004 


8.078 


8, 


.218 


0, 


.482 


0.306 


0.330 


7, 


.13 


it 


0.02 


6, 


.85 


it 


0.02 


6, 


.83 


it 


0.02 


-1.49 


0.37 


i 
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Table 1. Continued 



name 


HIP 




(mas) 




(K) 


e±Ae 

(mas) 


msol 




V 


B-V 


V-Rc 


{R-I)c 






J 








H 






Ks 


[Fe/H] 


[a/Fe] 




HD 


97320 


54641 


17. 


.77 


± 





.lb 


6142 ± 51 


0.180 


± 


0.003 


8.042 


8. 


.161 


0.481 


0.295 


0.312 


i-r 

7, 


1 A 

.14 


1 

± 


0.02 


6. 


.o7 


1 

± 


0.03 


6. 


.79 


± 


0.02 


-1.28 


0.30 


i 


HD 


98281 


55210 


45. 


.48 


± 


1. 


fin 
.UU 


5416 ± 58 


0.353 


± 


0.008 


7.125 


7. 


.275 


0.747 


U.4io 


0.375 


5, 


.yi 


1 

± 


0.02 


5. 


.57 


1 

zh 


0.03 


5. 


.46 


± 


0.03 


-U.2U 


n C\A 
U.U4 


b 


HD 


103072 


578DO 


o9. 


.48 


± 


1. 


.58 


5058 ± 52 


0.255 


± 


0.006 


8.129 


8 


.400 


0.850 


0.501 


U.4zD 


6 


.7o 


1 

± 


0.03 


D. 


OK 

.oD 


1 

± 


0.02 





.27 


it 


0.02 


n Qn 
-U.oU 


0.08 


( 

I 


HD 


105671 


59296 


49. 


,72 


± 


1, 


.00 


4604 ± 53 


0.339 


± 


0.008 


7.918 


8 


.440 


1.135 


0.660 


0.545 


6 


A A 

.44 


1 

± 


0.02 


5. 


.91 


1 

± 


0.04 


5 


.76 


± 


0.02 


0.20 


0.00 


c 


HD 


106156 


59572 


32. 


,30 


± 


1, 


.02 


5428 ± 45 


0.259 




0.005 


7.788 


7 


.918 


0.800 


0.426 


0.365 


6 


.57 


1 

± 


0.02 


6. 


.23 


1 

± 


0.04 


6 


.12 


ih 


0.03 


0.23 


-0.03 


b 


HD 


106516 


59750 


44. 


,34 


± 


1, 


.01 


6296 ± 135 


0.431 


± 


0.019 


6.038 


6 


. 104 


0.466 


0.286 


0.285 


5 


.13 


1 
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Table 1. Continued 
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6 


.77 


± 


0, 


,05 


6, 


.70 


± 


0, 


.02 


-0.62 


0.17 


n 


BD+51 1696 


57450 


13, 


.61 


± 1, 


.54 


5721 


it 


62 


0, 


.094 


± 


0.002 


9.760 


9.916 


0, 


.568 


0.354 


0.361 


8, 


.68 


± 


0, 


.02 


8, 


.38 


it 


0, 


.03 


8, 


.31 


it 





.02 


-1.48 


0.33 


P 


BD-01 306 


10449 


16, 


.17 


± 1, 


.34 


5771 


± 


50 


0, 


.133 


± 


0.002 


8.969 


9.084 


0, 


.583 


0.341 


0.355 


7, 


.90 


± 


0, 


.02 


7, 


.59 


it 


0, 


.04 


7, 


.52 


it 


0, 


.03 


-0.96 


0.36 


b 


BD-11 4126 


80295 


20 


.67 


± 2, 


.17 


4788 


± 


36 


0, 


.121 


± 


0.002 


9.986 


10.383 


1, 


.023 


0.585 


0.482 


8, 


.58 


± 


0, 


.03 


8, 


.06 


it 


0, 


.03 


7, 


.95 


it 





.03 


0.20 


-0.01 


q 


CD-35 360 


5004 


16, 


.28 


± 1, 


.76 


5227 


± 


32 


0, 


.099 


± 


0.002 


10.040 


10.266 


0, 


.758 


0.442 


0.435 


8, 


.82 


± 


0, 


.02 


8, 


.36 


it 


0, 


.03 


8, 


.28 


it 


0, 


.02 


-1.15 


0.39 


i 


CD-41 6367 


54538 


16, 


.39 


± 1, 


.24 


5219 


± 


45 


0, 


.125 


± 


0.002 


9.540 


9.738 


0, 


.785 


0.437 


0.408 


8, 


.29 


± 


0, 


.03 


7, 


.88 


it 


0, 


.04 


7, 


.76 


it 


0, 


.03 


-0.32 


0.09 


e 



Eflfective temperatures (T^ff) and angular diameters (0) are those computed via IRFM as described in Section 4. Apparent bolometric magnitudes {mBol) are obtained as described 
in Section 7, where the absolute bolometric magnitude of the Sun Mg^i^Q = 4.74. Source of metallicities : (a) Mishenina et al. 2002, 2004; (b) Valenti & Fischer (2005); (c) Pompeia, 

Barbuy &: Grcnon 2002, 2003; (d) Santos ct al. (2005); (e) Favata, Micc'la & Sciortino (1997); (f) Gray ct al. (2003); (g) Paulson, Sncdcn & Cochran (2003); (h) Santos ct al. (2004) 
(Mg abundances from Beirao, Santos et al. 2005); (i) Gratton et al. (2003); (j) Thoren & Feltzing (2000); (h) Clementini et al. (1999); (k) Fulbright (2000); (1) Brewer & Carney (2006); 
(m) Luck & Heiter (2005); (n) Tomkin & Lambert (1999); (o) Yong, Lambert & Ivans (2003); (p) Bai et al. (2004); (q) Castro et al. (1997). 



